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ABSTRACT 



We establish an isoperimetric inequality for the fundamental tone (first nonzero eigenvalue) of 
the free plate of a given area, proving the ball is maximal. Given r > 0, the free plate eigenvalues 
co and eigenfunctions u are determined by the equation AAu — tAu — uju together with certain 
natural boundary conditions. The boundary conditions are complicated but arise naturally from 
the plate Rayleigh quotient, which contains a Hessian squared term |D 2 it| 2 . 

We adapt Weinberger's method from the corresponding free membrane problem, taking the 
fundamental modes of the unit ball as trial functions. These solutions are a linear combination 
of Bcsscl and modified Bcsscl functions. 
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CHAPTER 1 

Plates and the isoperimetric problem 



Isoperimetric problems are about minimizing or maximizing a quantity subject to constraints. 
The classical isoperimetric inequality states that of all planar regions of the same perimeter, the 
disk has maximal area. Equivalently, of all regions of the same area, the disk minimizes perime- 
ter. The three-dimensional version can be observed physically quite readily - many mammals 
sleep curled up in a ball, while keeping their volume the same, to minimize surface area and 
hence heat loss. 

Many physical quantities satisfy isoperimetric-type inequalities. The goal of this thesis is to 
prove an isoperimetric result for a free plate under tension with unconstrained edges: of all such 
plates having the same area, the disk has the highest fundamental pitch. 

Researchers have investigated and proved isoperimetric inequalities regarding frequencies 
of vibration in related situations. Lord Rayleigh conjectured, and Faber and Krahn proved, 
that of all membranes of the same area with constrained edges, a circular drum produces the 
lowest pitch. Kornhauser and Stakgold conjectured the opposite bound for a membrane with 
unconstrained edges; this result was proven by Szego and Weinberger. This thesis generalizes 
their result to plates under tension. Plate problems are more difficult than membrane problems 
because they involve the bi-Laplacian rather than the Laplacian. 

Mathematical formulation 

We now develop the mathematical formulation of the free plate isoperimetric problem. Let SI 
be a smoothly bounded region in M d , d > 2, and fix a parameter r > 0. The "plate" Rayleigh 
quotient is 

r , L\D 2 u\ 2 +r\Du\ 2 dx , 
Q[u] = J "' f | | ' ' — ■ 1-1 

Here \D 2 u\ = (%2jk ^.jj 1 ' 2 is the Hilbert-Schmidt norm of the Hessian matrix D 2 u of u, and 
Du denotes the gradient vector. 
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Physically, when d = 2 the region D, is the shape of a homogeneous, isotropic plate. The 
parameter r represents the ratio of lateral tension to flexural rigidity of the plate; for brevity 
we refer to r as the tension parameter. Positive r corresponds to a plate under tension, while 
taking r negative would give us a plate under compression. The function u describes a transverse 
vibrational mode of the plate, and the Rayleigh quotient Q[u] gives the bending energy of the 
plate. 

From the Rayleigh quotient (1.1 1, we will derive in Chapter [2] the partial differential equation 



and boundary conditions governing the vibrational modes of a free plate. The critical points of 



(1.1) are the eigenstates for the plate satisfying the free boundary conditions and the critical 



values are the corresponding eigenvalues. The equation is: 

AAu- tAu = luu, (1.2) 

where uj is the eigenvalue, with the natural (i.e., unconstrained or "free") boundary conditions 
on dil: 

d 2 u 

Mu:= — =0 (1.3) 

Vu := r~ - div ao (P dn [{D*u)n}) - = (1.4) 

Here n is the outward unit normal to the boundary and divgji and grad^Q are the surface 
divergence and gradient. 



The eigenvalue equation (1.2) can also be obtained by separating the plate wave equation 

<j> tt = ~AA(f> + tA</>, 



by the separation </>(x,t) — u(x) cos(v / wt). The eigenvalue lu is therefore the square of the 
frequency of vibration of the plate. The quantities appearing as boundary conditions have 
physical significance as well. The expression Mu is the bending moment. As the plate bends, 
one side compresses while the other expands, leading to a restoring moment which must vanish 
at an unconstrained edge. 
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The problem 



We will prove in Chapter[2]that the spectrum of the free plate under tension is discrete, consisting 



entirely of eigenvalues with finite multiplicity: 



= ujq < u>i < uji < ■ ■ ■ — > oo. 



We also have a complete L 2 -orthonormal set of eigenfunctions u = const, u\, 1*2, and so forth. 

We call U\ the fundamental mode and the eigenvalue u)\ the fundamental tone; the latter 
can be expressed using the Rayleigh-Ritz variational formula: 



In general, the fcth eigenvalue is the minimum of Q[u] over the space of all functions u L 2 - 
orthogonal to the eigenfunctions uq, u\,. . . , ttfc-i. Because Uq is the constant function, the 
condition u J- uq can be written J^udx = 0. 

Let ST denote the ball with the same volume as f2. The main goal of this thesis is to prove 
the following theorem. 

Theorem 1. For all smoothly bounded regions of a fixed volume, the fundamental tone of the 
free plate with a given positive tension is maximal for a ball. That is, if r > then 



In the limiting case r = 0, the first d + 1 eigenvalues of fl are trivial because Q[u] =0 for 
all linear functions u. Thus we need the tension parameter r to be positive to get a nontrivial 
conjecture. 

The remainder of this chapter consists of a summary of the dissertation and then a brief 
history of related problems. 

We examine the behavior of the spectrum in Chapter[2j In particular, we prove the spectrum 
is comprised only of eigenvalues of finite multiplicity, with an associated complete eigenbasis. 
We obtain some bounds on the fundamental tone as a function of tension r in Chapter [3j where 
we also examine the fundamental tone in the extreme cases of infinite tension and zero tension. 
In addition, we derive the natural boundary conditions from the Rayleigh quotient. 

A discussion of ultraspherical Bessel functions appears in Chapter |4j along with a collection 




with equality if and only if Q is a ball. 



(1.5) 
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Figure 1.1: The fundamental tone Wi(fi*) of the disk (middle curve) and of another region wi(f2) 
(bottom curve). The dashed straight line is tangent to CJi(0*) at r = 0. 

of recurrence relations and facts about them that we will need to prove our main theorem. We 
put these facts to immediate use in Chapter [5j where we find the general form of eigcnfunctions 
for the ball and establish the angular dependence of its fundamental mode, in Theorem [2j 
Chapter [6] presents the proof of our main result, Theorem [I] 

In Chapter [7J we discuss the one-dimensional analogue of the free plate problem. Although 
there is no isoperimetric inequality in one dimension, we can prove an analogue of Theorem [2j 
that the fundamental mode of the free rod is an odd function about its midpoint. Also discussed 
are the eigenfunctions of the free rod under compression (r < 0). Building on that approach, 
In Chapter [8j possible future directions and generalizations of the main problem are discussed. 
Finally, in the Appendix, we gather some calculus facts used in prior chapters. 

Brief history of isoperimetric problems 

Isoperimetric problems for eigenvalues of the Laplacian have fascinated researchers for quite 
some time [3J [51 QUI [T71 [301 [531 \M\- in 1877, Lord Rayleigh [23J conjectured an isoperimetric 
inequality for the first eigenvalue of a fixed membrane: 

Ai(O) > Ai(fT) with equality if and only if fi is a ball. 

Here Xj is the jth eigenvalue satisfying the membrane equation — Au — Xu together with the 
boundary condition u = on dft. The above inequality was later proven by Faber [15] and 
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Krahn [35] and now bears their names. Kornhauser and Stakgold [53] conjectured in 1952 
the opposite result for the free membrane problem, —Au = fiu with du/dn = on dfl. In this 
case, the eigenvalues are = /Xo < /xi < (J>2 ■ ■ ■ ■ The lowest eigenvalue [io = corresponds to a 
constant eigenfunction. This mode does not vibrate; thus we call the next eigenvalue, n\, the 
fundamental tone. The Kornhauser-Stakgold conjecture sought to maximize the fundamental 
tone of the free membrane: 

Mi(^) ^ with equality if and only if f2 is a ball. 

This was proven by Szego (46j [48] and Weinberger [52] ■ Szego's proof uses conformal mapping 
and is only valid in two dimensions for simply connected regions. Weinberger's approach works 
for arbitrary domains in all dimensions. Our proof for Thcorcm[l]is inspired by Weinberger's ap- 
proach, namely using trial functions and demonstrating monotonicity properties of the resulting 
quotient of integrals. 

The bi-Laplacian operator appears in plate problems just as the Laplacian underpins mem- 
brane problems. The theory of the bi-Laplacian is not nearly so well developed as that of the 
Laplacian. For example, solvability of the biharmonic equation in Lipschitz domains with Neu- 
mann boundary conditions was established only a few years ago [51] . Furthermore, the maximum 
principle fails for the bi-Laplacian, even in one dimension. Even so, isoperimetric problems for 
plate eigenvalues have long been under consideration. Lord Rayleigh [43] conjectured an isoperi- 
metric inequality for the clamped plate with zero tension, A 2 u = au, where "clamped" refers 
to the boundary conditions u = du/dn = on <90. This isoperimetric inequality (like that of 
the fixed membrane) gives a lower bound on the first eigenvalue: 

ai(Q) > ax(n*). 

Work towards this result was begun by Szego [46] and continued by Talenti [49] (see also Mohr 
[31]). Nadirashvili [32], [33] proved the conjecture in two dimensions based on Talenti's work on 
rearrangement of elliptic partial differential equations. Ashbaugh and Benguria later extended 
the proof to three dimensions [8] . The problem remains open for dimensions four and higher, 
with a partial result by Ashbaugh and Laugesen [10] . and is open in all dimensions for r > 0. 
For an overview of work on the clamped plate problem, see [T71 Chapter 11, p. 169-174] and 
[22 p. 105-116] 
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There are other boundary conditions for the plate besides the natural and clamped conditions 
discussed so far. The simply supported plate is governed by A 2 — tAu = /3u with the requirement 
that u = on dQ; from here Mu = arises as a natural boundary condition. It is natural to 
conjecture an isoperimetric inequality for this problem too, but no work seems to have been 
done on it. 

There is also a body of work on plate problems that does not focus on isoperimetric in- 
equalities. General vibrating plate eigenvalue problems are discussed with experimental data 
by Leissa in |27) . including approximate solutions for the free rectangular plate. The buckling 
eigenvalues of a clamped plate have been considered by Payne [35J [33] ■ Supported plate work 
includes Payne 138. and Licar and Warner |28j . who examine domain dependence of plate eigen- 
values. The free plate without tension is considered in the same papers. Free plate work without 
tension also includes Nakata and Fujita, who establish upper and lower bounds on free plate 
eigenvalues in 34] . Payne considered the buckling problem for the clamped plate ([36, 38J, and 
in conjunction with Weinberger [32]). Kawohl, Levine, and Velte [21] considered the clamped 
plate under tension and compression; they viewed the eigenvalues as functions of the tension or 
compression parameter and established upper and lower bounds in terms of these parameters. 
Analogous results for free plates are in Chapter [3] of this dissertation. 

Isoperimetric inequalities for the Laplacian and related problems were first considered more 
than 130 years ago, and research in this field continues today. A number of problems regarding 
eigenvalues of the bi-Laplacian remain open, and I hope this dissertation helps provide insight 
for future research on plate problems. 
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CHAPTER 2 
The spectrum 



Our first task is to investigate the spectrum of the fourth-order operator associated with our 



Rayleigh quotient Q in (1.1). In this chapter we show there is only discrete spectrum, with an 
associated weak eigenbasis. We will then establish regularity of the eigenfunctions up to the 
boundary and derive the natural boundary conditions. 

In this chapter we will allow r to be any real number. We require fi C M. d to be smoothly 
bounded. 



The existence of the spectrum 

We consider the sesquilinear form 

p d 

a(u, v) = / V] u XiXj v XiXj + t(Du ■ Dv) dx 

in L 2 (fl) with form domain H 2 (fl). Note the plate Rayleigh quotient Q can be written in terms 
of a, with Q[u] — a{u,u) /\\u\\ 2 L2 . 

Proposition 1. The spectrum of the operator A associated with the form a(-, •) above consists 
entirely of isolated eigenvalues of finite multiplicity ujq < uj\ < UJ2 < ■ ■ ■ — > oo . There exists an 
associated set of real-valued weak eigenfunctions which is an orthonormal basis for L 2 (Q). 

Proof. By Cauchy-Schwarz, the form a(-,-) is bounded, and therefore continuous, on H 2 (Q). 
We will show the quadratic form a(u, u) is coercive; that is, for some positive constants c\ and 
C2, we have a(u, u) + ci|ju|j 2 > C2||u|| ^2^. By the boundedness of a on H 2 , this is equivalent 
to showing that the norm 

\\ u \\l = a(u,u) +c 1 ||u|| 2 , 
is equivalent to || • 11^2(5^), and hence a is a closed quadratic form on H 2 (Q). 
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Once we have coercivity and show that H 2 (£l) is compactly embedded in L 2 (fi), we conclude 
by a standard result (see e.g., Corollary 7.D p. 78]) that the form a has a set of weak 
eigenfunctions which is an orthonormal basis for L 2 (f2), and the corresponding eigenvalues are 
of finite multiplicity and satisfy 

w < w l < ' ' ' < w n 00 as 71 — ► 00- (2-1) 

We first show that H 2 (Q) is compactly embedded in L 2 (£l). Let -Br be a ball of radius R 
centered such that fl C Br; thus we have H 2 (Q) C L 2 (£l) C L 2 {Br) and can identify L 2 (f2) 
with a closed subspace of L 2 {Br). Because fl is smoothly bounded, we can extend all functions 
in H 2 (yt) to functions in H 2 (B2r). The extension map is linear and bounded, so H 2 (Vt) can 
be identified with a closed subset of H 2 (B2r). The space H 2 (B2r) is compactly embedded in 
L 2 (Br) by the Rellich-Kondrachov Theorem; thus because the space H 2 (Vt) is a closed subset 
of H 2 (B2r), we have that H 2 {yi) also compactly embedded in L 2 (Br). Since L 2 (fl) is a closed 
subspace of L 2 (Br) containing H 2 (il), we have that H 2 (fl) is a compact subspace of L 2 (fl). 
(See, e.g., [2] for extension theorems and the Rellich-Kondrachov Theorem.) 

We next show coercivity of the form a. For r > 0, coercivity is easily proved: 

a(u,u) +t\\u\\ 2 > \\D 2 u\\ 2 + t\\Du\\ 2 + t\\u\\ 2 
> min(r, l)||u||#2, 

where all unlabeled norms are L 2 norms on f2. 

To prove coercivity when r < 0, we must somehow arrive at a positive constant in front of 
the \Du\ 2 term. We cannot use Poincare's inequality on the \D 2 u\ term as this will introduce 
terms involving the average value of Du. Instead, we will exploit an interpolation inequality. 

By Theorem 7.28 of [THl p. 173], we have that for any index 1 < j < n and any e > 0, 

HQ^Hi^n) < elMla^n) + C^IMI^n) (2-2) 
with C = C(O) a constant. Replacing e by e/d and summing over j, we see 

\\D 2 u\\h > Q - 1) - (§ + l) II*- 
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Fix S £ (0, 1). Let K > 0. Then 




+ IN 



We can choose our e small and our if large so that the minimum is positive, which proves 
coercivity. For example, for S = 1/2, we need to take 



We now have that the form a is coercive for all r G E. Now suppose u is a weak eigenfunction 
corresponding to eigenvalue uj. Because u> and r are real- valued, by taking the complex conjugate 
of the weak eigenvalue equation we see that u is also a weak eigenfunction with the same 
eigenvalue. Thus the real and imaginary parts of u are both eigenfunctions associated with u>, 



Note that for any bounded region Q and all real values of r, the constant function solves the 
weak eigenvalue equation with eigenvalue zero. For all nonnegative values of t, the Rayleigh 
quotient is nonnegative for all functions and so = wq < w± < ■ ■ ■ . When r = 0, the coordinate 
functions x\, . . . ,Xd are also solutions with eigenvalue zero, and so the lowest eigenvalue is at 
least d + 1-fold degenerate, as noted in the introduction. Taking instead r > 0, the Raleigh 
quotient shows that the fundamental tone ui\ is positive, and so we have: 



Regularity 

We aim to establish regularity of the weak eigenfunctions by appealing to interior and boundary 
regularity theory for elliptic operators. 

Proposition 2. For any r € R and smoothly bounded £1, the weak eigenfunctions of the operator 
A are smooth on f2. 

Proof. Let u be a weak eigenfunction of A with associated eigenvalue lu; by Proposition [l] we 
have u € T>(a) — H 2 (fl). Then by a theorem in [33 p 668], we have u G H k (il) for every 




and we may choose our eigenfunctions to be real-valued. 



□ 



= ljq < uj\ < ui<x < • • • < u> n — »■ oo as n oo. 
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positive integer k. Thus we have u £ H k (fl) for all k € Z + , and so u € C°°(f2). 

Regularity on the boundary follows from global interior regularity and the Trace Theorem 
(see, for example, [50, Prop 4.3, p. 286 and Prop 4.5, p. 287.]). Thus we have u £ C°°(fi), as 
desired. □ 



The Natural Boundary Conditions 

In this section, our goal is to derive the form of the natural boundary conditions necessarily 
satisfied by all weak eigenf unctions. 

In the case of the free membrane, the weak eigenfunctions u are smooth on VL and satisfy 

Du ■ Dip dx — [iu(j> dx — 
n Jn 

for all <p £ iJ 1 (SI) . By the Divergence theorem, we obtain 

/ 4>^dS- [ (j)(Au + fxu) dx = 0. (2.3) 
Jon on J n 

Because 4> € -ff 1 (£1) is arbitrary, we may consider those </> with compact support in f2. Then 
the boundary terms vanish, and in order for the remaining integral to be zero for all such <fi, we 



must have Au + fiu = almost everywhere on f2. Hence equation (2.3) says 

Jan on 

In order for this surface integral to be zero we find must have du/dn — on <9f2. Thus u satisfies 
the eigenvalue equation — Au = uu on SI and the Neumann boundary condition 

5" = o on an. 

on 

We will use the same approach to derive the natural boundary conditions for the free plate. 
The natural boundary conditions are rather complicated in higher dimensions, and so we state 
the two-dimensional case first. The boundary conditions in this case have been known for some 
time: see, for example, [53] 

Proposition 3. (Two dimensions) For fl C M 2 , the natural boundary conditions for eigenfunc- 
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Hons of the free plate under tension have the form 

y du d(Au) d / d 2 u ^du 
dn dn ds \dsdn ds 

where n denotes the outward unit normal derivative, s the arclength, and K the curvature of 

dn. 

We also look at one example of the natural boundary conditions for a region with corners. 
Notice that an additional condition arises at the corners! 

Proposition 4. (Rectangular region in two dimensions) When fl C M. 2 is a rectangular region 
with edges parallel to the coordinate axes, the natural boundary conditions for eigenf unctions of 
the free plate under tension have the form 



a u , , 

——, - = at each edge 
on 

du d 3 u d(Au) 

T x » — = on each ed 9 e 

dn ds z dn dn 

u xy = at each corner 



where n and s indicate the normal and tangent directions. 

Finally, we state the natural boundary conditions for a smoothly-bounded region in higher 
dimensions: 

Proposition 5. ( General) For any smoothly bounded Q, the natural boundary conditions for 
eigenfunctions of the free plate under tension have the form 

d 2 u 

Mu := — ^ = on dfl, 

dn 2 - 

du ,. r ,„ s i \ dAu 

Vu := t- 



dn 



Pan [(D 2 u)n\ J - — =0 on dQ, 



where n denotes the normal derivative and divg^ is the surface divergence. The projection Pqq 
projects a vector v at a point x on dfl into the tangent space of dQ at x. 

Proof of Proposition^ Our eigenfunctions u are smooth on Q, by Proposition 2 and satisfy the 
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weak eigenvalue equation a(u, <j>) — uj(u, </>)l 2 (si) = for all <fi e H 2 (fl). That is, 
/ ^ u XiX . 4> XiXj + tD4> ■ Du - u)U(j) I dx = 0. 

As in the membrane case, we make much use of integration by parts. Let n denote the 
outward unit normal to the surface dfl. To simplify our calculations, we consider each term 
separately. 

The gradient term only needs one use of integration by parts: 

/ Du-Dcj)dx= f 4>^dS- f <f>(Au)dx. 
Jn Jan v n Jn 

The Hessian term becomes: 

= J (d<P- ([D 2 u)n^j - <t>^j^j dS + J (A 2 u)cf> dx, 

after integrating by parts twice. 

We wish to transform the term involving Dcj) in the above surface integral using integration 
by parts. Because we are on dfl, we must treat the normal and tangential components separately. 
We can then use the Divergence theorem for integration on dfl. 

We note that the surface gradient grad an equals D — nd n when applied to a function (like 
</>) that is defined on a neighborhood of the boundary. Thus grad ao </> gives the tangential part 
of the Euclidean gradient vector. Hence, 

/ D4> ■ ((D 2 u)n) dS 
Jan ^ ' 

= L ("^ + ■ ("£ + P »" [< D2 »)"]) " S 

by the Divergence Theorem on the surface dfl. Here (•, -)gn denotes the inner product on the 
tangent space to dfl. Recall Pqq projects a vector at a point x on dfl onto the tangent space of 
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dfl at x. 

Thus for u an eigenfunction associated with eigenvalue w, we see 

= J <j)(^A 2 u — tAu — uju^i dx 

f dcf)d 2 u , / du dAu t r 2 , \\ 

As in the membrane case, this identity must hold for all <f> 6 H 2 (Q). If we take any compactly 
supported 0, then the volume integral must vanish; because <f> is arbitrary, we must therefore 
have A 2 u — tAu — urn = everywhere. Similarly, the terms multiplied by <fi and d(f>/dn must 



vanish on the boundary. Collecting these results, we obtain the eigenvalue equation (1.2) and 



natural boundary conditions of Proposition [5j □ 

Proof of Proposition^ Here d — 2; take rectangular coordinates (x,y). We parametrize dfl 
by arclength s and define coordinates (n,s), with n the normal distance from dCl, taken to be 
positive outside fi. Write n(s) and t(s) for the outward unit normal and unit tangent vectors 
to the boundary. Then Pqq \f\n + /2^] = fit and the operators div^Q and grad^Q both simply 
take the derivative with respect to arclength s. That is, for a scalar function /(s), and taking 
t(s) to be the tangent vector to the surface, we have 

grad ao /( S ) = f(s) and div dn (f( S )t(s)) = f'(s). 



and so we may write 

div^ (^n [(D 2 u)n] ) = ^-t T (D 2 u)n. 

The tangent line to d£l at the point (0, s) in our new coordinates forms an angle a = a(s) 
with the x-axis (see [53j p. 230]); the curvature of d£l is given by K(s) = a'(s). Then in 
rectangular coordinates, the unit tangent vector is (cos a, sin a), and the outward unit normal 
is (sin a, — cos a). Thus we have 

—t T (D 2 u)s — d s ( sin a cos a(u xx — u yy ) + (sin 2 a — cos 2 a)u xy \ . 
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By [531 P- 233], on dfl under our change of coordinates, we have 

u X x = Unn sin 2 a + u ss cos 2 (x + 2u ns sin a cos a + Ku n cos 2 a — 2Ku s sin a cos a 
w TO = u nn cos 2 a + m ss sin 2 a — 2u ns sin a cos a 4- -fCu„ sin 2 a + 2Ku s sin a cos a 
= —u nn cos a sin a + u ss cos a sin a + w ns (sin 2 a — cos 2 a) 
+ if u n cos a sin a — Ku s (sin 2 a — cos 2 a) . 

So after simplification, 

This together with the results of Proposition [5j yields the form of Vu given in Proposition [3j 
Mu is unchanged, and so this completes the proof. □ 

Proof of Proposition^ Our previous findings do not completely apply because dft has corners, 
although our argument proceeds similarly. For convenience of notation, we will take fi to be 
the square [0, l] 2 . 

The Hessian term gives us a condition at the corners. In particular, after integrating by 
parts twice, we have: 



Since 



and 
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we obtain 



,dA 



^(^XXXX H~ 2>U XX yy + Uyyyy^ 



dA 



+ y {^ Uj yy ( t ) y ~ ( \ > {'^ u xxy + u yyy)^J 



x=l 
x=0 

y=0 



+ 2u 



■'•y 



a;=0 



y=0 



Because the Divergence Theorem does apply to regions with piecewise-smooth boundaries, 
the gradient term is the same as in the smooth-boundary case. The final term above is the only 
term that depends only on the behavior of u and <f> at the corners; arguing as before, we obtain 
the eigenvalue equation and natural boundary conditions, with the additional condition 



= u 



xy'-: 



1 

x=0 



1=0 



That is, we must have u xy — at the corners. 



□ 



Example: natural boundary conditions on the ball 

When is a ball, we can simplify the general boundary conditions. 

Proposition 6. (Ball) The natural boundary conditions in the case 
radius R, are 



, the ball of 



Mu := u rr = 
Vu 



:= ru r - ^A s (u r - ^) - (Au)r = 



atr = R, (2.4) 
atr = R. (2.5) 



Proof. When fi is a ball, the normal vector to the surface at a point x is n = x/R. Then the 
ith component of (D 2 u)n is given by 



EXj 
u XiXj — 

3 = 1 
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and can be rewritten as 

d 



1 d_ 
R dxt 



Xj-U 



Therefore, 

(D 2 u)n = D ( Du ■ - - ^ 

K ' V R R/ 

Then the projection Pqq takes the tangential component of the above gradient vector, and so 

Pan [{D 2 u)n\ = grad 0Bd(fl) (Du • -| - . 

We know divgngrad,^ = AgQ by definition. For the ball of radius R, we have A 9B d^ — ^jAg. 
The operator Ag is the spherical Laplacian, consisting of the angular part of the Laplacian. It 
satisfies the identity A = Jp + + ^ A s . 

Thus 



divonPdn [{D u)n\ = A dB d {R) yDu ■ — - 



by noting that • n = du/dr. The boundary conditions of Proposition [5] then simplify to ( 2.4 1 



and (2.5), as desired. □ 
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CHAPTER 3 



The fundamental tone as a function of tension 



Fix the smoothly bounded domain Q. We will estimate how the fundamental tone ui\ = wi(t) 
depends on the tension parameter r, for use in the proof of Theorem [l] We will also study in 
this chapter the behavior of ui% in the extreme cases as r — > and r — > oo. 



First we note that the Rayleigh quotient (1.1 1 is linear and increasing as a function of r 



Our eigenvalue wi(r) is the infimum of Q[u] over u £ H 2 (fl) with J^udx — 0, and thus wi(r) 
is itself a concave, increasing function of r. 

Next, we will prove uj\{t)/t is bounded above and below for all r > 0. Recall fi% is the 
fundamental tone of the free membrane. 

Lemma 3.1. For all r > we have 

tih<Ui(t )<r - r | 1 , (3.1) 

J n \x-x\ 2 dx 

where x = J a xdx/\Cl\ is the center of mass of VI. In particular, when ft is the unit ball, 

TfJ,i < wi(r) < r{d + 2). (3.2) 



Furthermore, the upper bounds in (3.1 1 and (3.2 I hold for all r G 



These bounds are illustrated in Figure [XT] 

Proof. To establish the upper bound, take the coordinate functions as trial functions: Uk = 
Xk — Xk, for k = 1, . . . , d. Note J n Uk dx = by definition of center of mass, so the Uk are valid 
trial functions. All second derivatives of the Uk are zero, so we have 

u 1 {T)<Q[u k ] = $ * r{DUk?dX J - ldX 



In u l dx J n (x k - x k ) 2 dx' 
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Clearing the denominator and summing over all indices k, we obtain 

Wi(r) / \x — x\ 2 dx < r\il\d, 
Jn 

which is the desired upper bound. When il is the unit ball, note J n \x\ 2 dx = \Vt\d/ (d + 2). 
Now we treat the lower bound. Let u £ H 2 (£l) with J Q udx = 0. Then 

r , t L \Du\ 2 dx 

Q[u] > J "' ' > r Ml 

in u dx 

by the variational characterization of fi\. Taking the infimum over all trial functions u for the 
plate yields wi(r) > Tfi\. □ 

Note that Payne |37j proved linear bounds for eigenvalues of the clamped plate under tension. 
Kawohl, Levine, and Velte [3T] investigated the sums of the first d eigenvalues as functions of 
parameters for the clamped plate under tension and compression. 

We can also prove another linear upper bound on uj\, which is just a constant plus the lower 
bound in Lemma 13. II 



ui < <7(n)+rMi, 



Lemma 3.2. For all r£l, 



where the value 



, . L\D 2 v\ 2 dx 
Jn v dx 
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is given explicitly in terms of the fundamental mode v of the free membrane on O. 

Proof. Let v be a fundamental mode of the membrane with Av = —fi\v and J^vdx = 0; the 
membrane boundary condition is du/dn = on dQ. Then by the variational characterization 
of eigenvalues, 

w x (t) < Q[v] = C(Sl) + tQ m [v] - C(fi) + Tin, 
as desired. □ 



tUj (r) 




r 



Figure 3.2: The fundamental tone of the disk (solid curve) together with the upper bound of 
Lemma 3.2 (top dashed line) and the lower bound of Lemma 3.1 (bottom dashed line). 



Infinite tension limit 

A plate behaves like a membrane as the flexural rigidity tends to zero, that is, as r = (tcnsion/flexural rigidity) 
tends to infinity. For the fundamental tone, that means: 

Corollary 3.3. For the fundamental tone of the free plate, 



— > fi\ as t — y co. 



Proof. By Lemmas 3.1 and 3.2 we have 



r r 



Let t — v oo. 



□ 
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The eigenfunctions should converge as r — > oo to the eigenfunctions of the free membrane 
problem. Proving this for all eigenfunctions seems to require a singular perturbation approach, 
which has been carried out for the clamped plate in |14) , but we will not need any such facts for 
our work. For the convergence of the fundamental tone of the clamped plate to the first fixed 
membrane eigenvalue, see |21j . 

Vanishing tension limit; moment of inertia 

At r = 0, the lowest eigenvalue is zero and has multiplicity d + 1 since Q[u] — for any linear 
function u. We will establish a relationship between the (scalar) moment of inertia In of our 
region Q and the derivatives at r = of the first d nontrivial eigenvalues wi(r), . . . , w^t). 

Lemma 3.4. For t > 0, we have 



We will not need this result later, except as motivation for some conjectures. 
If uij (r) is differentiable from the right at t = then we deduce the following bound involving 
the derivatives of eigenvalues with respect to r: 



Proof. We assume our plate has its center of mass at the origin, so that the scalar moment of 
inertia Iq may be expressed as 



where Mq — J„ xx T dx is the moment matrix whose entries are given by = J n XiXj dx. 

Take a(., .) to be the sesquilinear form from Chapter 2. As in [H] p. 99], we define the 
inverse trace tr _1 [L m ] of the m-dimensional space 
L m c{»€ H 2 {VL) : f n vdx = 0} by 



d 






or 





where the Wj form a basis of L m satisfying the orthonormality condition 
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a(wj,Wk) — Sjk- Then we have, again by the variational characterization 

^ — ^— = max { tT~ 1 \L d \ : L d C H 2 (Q) with / v dx = Vw £ L d \ . 



£ 

3= 



Considering the coordinate functions Xj, we see a(xj,Xk) = T\Cl\5jk and 



/n x 3 = ^' so we ma y take w i ~ x il V T l^l> an d our variational characterization gives us 

3=1 3V ' 3=1*'" 

The righthand side is simply the trace of Mq and hence is equal to our scalar moment of inertia 
la- □ 
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CHAPTER 4 



Ultraspherical Bessel functions 



We must examine properties of d-dimensional ultraspherical Bessel functions, for they provide 
the eigenfunctions on the unit ball for dimensions 2 and higher, in the next chapter. For more 
information on Bessel functions, see [TJ p. 358-389]. For more information on spherical and 
ultraspherical Bessel functions, see [U p. 437-455] (d = 3 only) and [30] (all d > 2). 



Definitions 

The Bessel function J u (z) is defined by the power series 

{-l) k (Z\^+v 



Jv{z) -^k\T( V + k + l) (2) 



o fc!r(^ + fc + i) 

and is hence analytic. This function solves the Bessel differential equation, z 2 w" + zw' + (z 2 — 
v 2 )w = 0. 

For higher dimensions, we need to consider spherical (d — 3) and ultraspherical (d > 4) 
Bessel functions ji(z), defined by: 

ji(z) = z~ s J s+l (z) 
d-2 

with s — . 



Such functions solve the equation 

7" + {d- l)zw' + (z 2 - 1(1 + d - 2)) w = 0. (4.1) 



z 2 w" 



Analogously, the modified Bessel function I v (z) is given by the power series 
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and solves the modified Bessel equation z 2 w" + zw' — (z 2 + v 2 )w = 0. We define the higher- 
dimensional analog ii(z) as follows: 

k{z) = z~ s I s+ i(z) 
d-2 

with s = — - — . 



Such functions solve the equation 

z 2 w" + (d- l)zw' - (z 2 + l{l + d-2)^w = 0. (4.2) 

Recurrence Relations and power series 

The Bessel functions Jy and l v have a number of useful recurrence relations. Those listed below 
are taken from [H P- 361, 376]. 

— J v (z) = Jv-l(z) + J v+ i(z) 

z 

J'»( z ) = -Jv{z) - J V +l{z) 

z 

V 

= Jv-\\Z) - -Jv(z) 

—I v (z) = h-i{z) - I„+i{z) 
I' v {z) = - z I v {z) + I v+l {z) 

= I v -i(z) - -I v {z) 

z 

From these we also have recurrence relations involving second derivatives: 

V 2 - v \ , , % 1 



= (- 



- 1 J v (z) + -J v+ i{z) 



z 2 ) ' z 



+l)I v {z)--J v+l {z) 



z 



The ultraspherical Bessel functions have similar recurrence relations, all of which follow from 
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the definition and application of the corresponding ordinary Bessel recurrence relations: 



d - 2 + 21 . , . . , , . , , 

3i{z) = 3i-\{z) + 3i+\{z) 

z 

fi(z) = -ji(z) -ji+i(z) 



. l+d-2 . 
: 3i-i{z) ji{z) 



d - 2 + 21 . 



k{z) = ii-i(z) - ii+i{z) 

i'l(z) = -il{z) + h+i(z) 
z 

= H-l{z) li{z) 

z 



(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 
(4.8) 



Note that if we take d = 2, each of these simplifies to the corresponding relation for Bessel 
functions. 

We also have recurrence relations for the second derivatives: 



fi'(z) 



l 2 -l A ., , d-l . , . 
l 2 -l A . , , d-l. , , 

— 1- 1 n{z) ii+i{z). 



(4.9) 
(4.10) 



Again, when d = 2 each recurrence relation simplifies to its two-dimensional analog. 

We may also write a power series for the ultraspherical Bessel functions ji(z) and ii(z) using 
the series for the corresponding J s+ i and / s +;: 



(4.11) 
(4.12) 



By examining the power series (4.12 1, it is immediate that ii(z) and its derivatives are all positive 



on (0, oo). Since the terms of the power series for j\ and ii are the same up to a sign, we also 
have that the derivatives of ji are dominated by those of i; : 



]\ m \z) <i\ m >(z) form>0, z > 0, 



■ (m) 



(4.13) 



with equality only at z = 0. 
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Other needed facts 

To prove our main result, we will need several facts about Bessel functions and their derivatives. 
We begin with a result on the zeroes of the j[{z). 

Proposition 7 (L. Lorch and P. Szego, [3D])- Let p^k denote the kth positive zero of j',(z). 
Then for d > 3 and I > 1 , 

l(d + 2l)(d + 2l + 2) 2 

< feil ) <i(d + 21). 

In particular, for p\ y \ the first zero of j[, we deduce 

d < p\ x < d + 2. 

T/izs inequality holds for all d > 2. 
Recall s = (ef-2)/2. 

Lemma 4.1. XTie functions ji and J s+ i have the same sign. In particular, for I = 1, . . . , 5 and 

anj/ d>2, we have ji(z) > /or z < pi\. 

Proof. The first statement is immediate from the definition of the ultraspherical Bessel functions. 
For the second statement, we appeal to established facts of Bessel functions. If we write for 
the first nontrivial zero of the Bessel function Ji(z). It is a well-known fact that Ji{z) is positive 
on (Ojj^i) and the zeroes are increasing in I for I > 1. Because Ji(z) — at z — and 
with no zeroes between, we have the same for ji(z) and thus the first root of j'i(z), px,x, 
lies between and Therefore for any d > 2 and any I > 1, we have Ji(z) > and hence 
ji(«)>0on(0 > pi, 1 ]. □ 

Lemma 4.2. We have j[ > on (0,pi.i). 

Proof. This follows from the observation that ji(z) > on (0, and the definition of pi i- □ 
Lemma 4.3. FFe /iawe j 2 > on (0,pi.i]- 

Proof. Let p2,i denote the first zero of j' 2 . By Proposition [7J p\ x < d + 2 and 

a 2(rf + 4)(d + 6) 

Then p| x -pf^ > (rf 2 + 8d + 28)/(10 + d) > 0, so that ,j 2 > on (0,pi,i]. □ 
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Lemma 4.4. We have j" < on (0,pi,i]. 



Proof. We see that 



il(*) = ^j2(z)-Jl(*) 



by (4.9) 



1 



h(z)-j' 2 (z) 



by (4.5) with I = 2. 



Since both j 2 , j' 2 are positive on (0,pi ; i] by the previous lemmas, we obtain j" on that same 
interval. □ 

Lemma 4.5. We have jj 4 ' > on (0,1*1,1]. 



Proof. We have by (4.9) that 



3x{z) = -jx{z) + d —^j 2 {z), 



and so 



d- 1 .„ 2(d- 1) ., 2(d- 1) 



L h(z) 



. . . 2(d- 1) . , , d 2 - 1 . . . 



(4.14) 



by (4.9) with Z = 1 and Z = 2, and (4.4) with Z = 2. When d = 2, this becomes 



# } W = (i- J)iiW + (^-| , ]./,u. 



(4.15) 



by (4.3) with I = 2. For any d, (4.14) gives us 



ti\z)~h{z) + {^-l)Hz) 



15 \ . . . d-4 . , . 



by (4.3) with Z = 3 



15(d + 6) 10\ . . . / 15 \ . 



by (4.3) with Z = 4. 



by (4.3) with / = 2 (4.16) 



(4.17) 

(4.18) 
(4.19) 

(4.20) 



When d — 2, then the first term of (4.16) is nonnegative on (0,pi,i] by Lemma 4.3 The 
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function j'3 is positive on (0,^1^] by Lemma 4.1 note that since d = 2, we have js(z) = Ja(z). 
Thus we have j[ (z) > when z G (0, n (0,pi,i]. However, p^i > y/3, so we have only 
established positivity on (0, y/S\. 



To establish positivity on (\/3,Pi,i] we turn to (4.15). The first term is certainly positive on 



(V^iPi,i]- The second term is positive when both J 2 > and z < \/6- Because ps 1.84 for 
d = 2, we have pi 1 < 6 and we are done. 



When d = 3 and d = 4, we again examine (4.16). Then Lemma 14. II together with the 



argument above give us > on (0, \/d 2 — 1) D (Ojjj^i]. By Proposition [7] we have pi.i < 
V d + 2, which for d = 3 and 4 is less than ydF=l, thus proving the lemma for these <i. 



For dimensions d > 5, we turn to (4.18). The second term is positive on (o,pi 1] for all d > 4 
by Lemma 4.1 Since pn < y/d + 2 and \/d + 2 < \/T5 for d < 13, we conclude j[ (z) > on 



(0,pi : i] for d = 5, . . . , 13. 

Finally, suppose d > 14 and z € (0,f>i,i]. If z 6 (0, \/T5], then j[ (z) > as above. If 



z > vl5, then we examine (4.20). Here the first term is nonnegative on [Vl5,Pi,i]. The 



non-Bessel factor of the second term is positive on (0, y §(d + 6)] and hence on (0,^1,1] 
Let dk denote the coefficients of the series expansion for i'{(z), so that 



□ 



i" (*) = 2j(-i) fc 4* and i"{z) = yy k z 2 

k=l 



fe=l 



by (4.11) and (4.12), where 



dfe = 



2k + 1 



(fc- l)!r(fc + l + d/2) 



2l-2fe-d/2 



Lemma 4.6. VFe have the following bounds: 



-diz + d 2 z 3 > j"(z) 
6 

d x z + -d 2 z 3 > 
5 



for all z e [0, v/3(d + 2)/(d + 5) 



/or 2; € 



0,V3 



Proof. Let 



d 



fe+X 



2fc + 3 



d fc 2/e(2fc + l)(2fc + d + 2)' 
It is easy to show that Ck is decreasing for k > 1. 
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We use the series expansion to first prove the following upper bound on f{(z) for z > 0: 

oo 

(-d lZ + d 2 z 3 ) - j»(z) = J2(-^) k+1 d k z 2k - 1 

k=S 
oo 

= ^(l-c fc z 2 )4^- 1 

k=S 
fcodd 

oo 

fc=3 
/codd 

since c k is decreasing in k. Hence (— d\z + d 2 z 3 ) — j'{{z) > when 

< z < 1/VcI = y/6(d + 4)/5, which is a larger range even than claimed in the first estimate 
in the lemma. 

For i"(z) we must take a slightly different approach. We will show that on [0, y/3], 



1 oo 

\d 2 z 3 >J2d k z 2k -\ 



fc=3 



and thus 



d lZ + %d 2 z z > i'[{z). (4.21) 
5 
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On [0, y/S\, note that 



^ k 2^(k-l)\T(k + l + d/2) \2J 



; _ 3 (fc-l)!r(fc+l + d/2) 



2fc + 1 /V3 



2fc-4 



< 



, (fe — l)!(fc + d/2)T(k + d/2) \ 2 

since z < v3 
2- d / 2 /z\3^ 2fc + l f^ k - 2 



T(3 + d/2) V2/ ^ (k- l)!(fc + d/2) \4 
since T(z) > 1 and is increasing on [2, oo), 

~ 5 f-^ fc! V4 

by the definition of c?2 and taking k t— > k + 1 

by the power series for e x 



5 

Thus we have obtained our desired bound on □ 
Bessel functions of the second kind 



Each of the Bessel equations (4.1 1 and (4.2 1 is a second-order differential equation, and so 
has another set of solutions. However, these functions are singular at the origin. We proved in 
Chapter [2] that the eigenfunctions are smooth; thus either these singular solutions do not appear 
in the eigenfunctions, or they appear in a linear combination such that the singular terms cancel. 
In Lemma |4.7[ we will prove that in fact there is no nontrivial linear combination that meets 
the smoothness condition. 



Ultraspherical Bessel functions of the second kind solve (4.1 ) and are defined by 

m (z) = z~ s N s+ i{z) 

with s — (d — 2)/2. Here N v (z) denotes a Bessel function of the second kind of order v. 
Each N v {z) is linearly independent of J„(z) (see, for example, [I] p. 358]), so ni(z) is linearly 
independent of ji(z). The functions Ni(z) are often written as Yi(z); we use Ni to avoid confusion 
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with the spherical harmonics Yi (9) . 

Ultraspherical modified Bessel functions of the second kind solve (4.2) and are defined by 



k l {z)=z- s K s+l {z) 
d-2 

with s 



2 ' 

where K v (z) denotes a modified Bessel function of the second kind of order v. As before, the 
ki(z) are linearly independent of the ii(z). 

We will need several properties of these functions. For orders v that are nonnegative integers, 
we have the following ascending series: (see, for example, [TJ p. 360, 375]) 

k=0 



1 oo 

ln(z/2) J„(z) - - V C v {k){-l) k z 2k+v , (4.22) 

7T * * 



2 

7T 7T 



fc=0 

m*) - \ £ (-i) fe (^/2) 2fe ^ 



1 oo 

+ \n{z/2)I v {z) + (-1)"- £ a(fc)(z/2) 2fe+ ^, (4.23) 



2 

A=0 



where the coefficients C v (k) are nonzero real number depending on z/ and fc as follows: 

V>(fc + 1) + ip{v + k + 1) 



C v {k) 



k\{v + k)\ 



Here - 0( n ) is the digamma function with "0(1) ~ —.577 and ip(n + 1) = ijj(n) + 1/n for n > 1. 

For positive noninteger orders ^, we have the following relations between Bessel functions of 
the first kind and second kind (see, for example, [U p. 358, 375]): 

= M^HW^M!) (424) 

sui(i/7r) 

, , 7T I_,.(z) — I,,(z) 

K v {z) = - " K \ 4.25 

2 sin(i/7r) 

When ^ is an integer, we have that J„ and J_„ are linearly dependant. 

Note that for all dimensions d > 2 and all nonnegative integers i, the functions n;(z) and 
fc;(z) are continuous for z > 0. 

In Chapter [3J we will find the exact solutions of our eigenvalue equation on the unit ball. 

30 



We will need the following lemma in order to show that Bcsscl and modified Bessel functions of 
the second kind do not appear in the radial parts of the smooth solutions. 

Lemma 4.7. Let a, b be positive constants, a < b. 

For d > 2 and all integers I > 2, there is no nontrivial linear combination 

R(z) = Ani(az) + Bki(bz) 

so that R(z) is smooth at z = 0. 

For d > 2 and 1 = 1, there is no nontrivial linear combination 

R(z) = Anx{az) + Bk\{bz) 

so that R(z) is smooth at z = with R(0) = 0. 

For d > 2 and I = 0, there is no nontrivial linear combination 

R(z) = An (az) + Bk (bz) 



so that R(z) is smooth at z = with R'(0) = 0. 

Proof. The Bessel functions n\ and ki are real-valued, so we may assume the constants A and 
B are real. Recall s — (d — 2)/2. We treat the cases of dimension d even and odd separately. 

[Part 1.] Let d > 3 be odd, so that s is an odd multiple of 1/2. Then for any nonnegative 
integer I, we have cos((s + 1)tt) = and sin((s + l)ir) = ±1. Thus by the definitions of ni and 



ki, and the identities (4.241 and (4.251, we have 



R(z) = Ani(az) + Bki{bz) 

= A{az)- s N s+l {az) + B{bz)- S K s+l {bz) 

= ± (-Alazyj-^az) + ~B(bz)~ s (l- s -i(bz) - I s+ i(bz))) 

y, ( z /2)2fc-2 8 -t 2 - a , n k _ 2s _ t _ k 2k - 2s -l\ 

^ k\T(k+l-s-l) V2 J 



k=0 



T -Bii(bz), 



by the power series expansions (4.11 ) and (4.12 1 for J v and l v . The terms in that contribute to 
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the singularity at the origin are: 



(z/2) 



2k— 2s— I 



^ kT(k + l- s-l) 



(i 



Bb lk 



1k-1s-l 



A{-l) k a 



k2k-2s-l 



So long as 2 — 2s — I = A — d — Z < 0, the function R(z) has at least two terms z 



-d+2-l 



and 



A-d-l 



that are singular at the origin (corresponding to k = and k = 1 in the above sum) 



Choosing A and -B so that the singular terms cancel, we must have 



7r fa\ d ~ 2 + l 



A 



2 \b 



d-4+l 



for k = 
for fc = 1. 



(4.26) 
(4.27) 



The constants a and 6 are positive; thus we must have A = B — for all d, I with 4 — tf — 1 < 0; 
this inequality holds for all odd d > 3 and nonnegative integers I except for the case of d = 3 
with I = or / = 1. 

When d — 3 and i = or I = 1, we only have one singular term (corresponding to k = 0). 



Assume B ^ 0, since if -B = then ^4 = by (4.26 1 and (4.27) and so the linear combination is 



trivial. Without loss of generality we may take B = 1. Then the linear combination 



R{z) 



|(|) ni(az) + ki(bz) 



is continuous at z = 0. Differentiating, we see 



fc=i 



(2fc-l-Z)(z/2) 2fc - 2 -' /tt 2fc _ x _, _ 7T 
2fc!r(fc + 1/2-0 V 2 2 



H-TgW^ftz). 



When / = 0, we have 



R'(0) = 



1 

2r(3/2) V2 



4r(3/2) 



2 v^y 



- |wj(o) 



noting io(0) = «i(0) = by < 4.7 ) and the series expansion (4.12) for ii(z). Because a ^ 6, we 
have i?'(0) 7^ for i = 0. Thus in order to have R{z) smooth at z = with R'(0) = 0, we must 
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take both A, B to be zero. 
When I = 1, we have 

Then since «i(0) = 0, we have 

R(0) = 



2r(i/2) V & 2 

which is nonzero. Thus in order to have R{z) smooth at z = with i?(0) = 0, we must take 
both A, B to be zero. 

[Part 2.] Ascending power series centered about zero with powers increasing by steps of 2 
and with lowest-order term Cz k will be represented by 0(z k ). Thus 4- 0(z k ) = 0(z k ^ 1 ) for 
k > 1 and 3jO(l) = O(z). 

Let d > 2 be even; then s = (d — 2)/2 is an integer. Thus by the definitions of ni and fc;, 



and the identities (4.22) and (4.23), we have 



R(z) = Ani(az) + Bk t (bz) 

= 2 _/^ 1 (g + /-fc-l)! /^2fe-2,-i /-A a2fe _ 2s _ ; + i? ( _ i)fc&2fc _ , 



fc! V2/ V 7T 2 

fc=0 v 

2/1 

+ ln(oz/2) — j^az) + (-l) s+ ' \n(bz/2)Bii(bz) + 0(z l ) (4.28) 

7T 

The singular contributions at z = come from the terms involving ln(z) when I — and when 
Z > from the sum 

*£' (■ + '-*-!)■ (£)»—' (=^-..-, + f , (4.29) 

fe=0 " ^ ^ ~ ' 

provided that at least one of s, I is nonzero. So long as 2 — 2s — Z = 4 — d — Z < and s 
and Z not both zero, the above sum contains at least two terms that are singular at the origin, 
corresponding to k — and k = 1. Choosing A and B so that both of these singular terms 
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cancel, we again find 



2 \bJ 

A= — B ( - 
2 \b 



d-2+l 



d-4+l 



for k = 0, 
for fc = 1. 



The constants a and 6 are positive; thus we must have A = B = 0. We have 4 — <i — I < for 
all even d > 6, for d — 4 when I > 1, and for d = 2 when / > 3. The remaining cases are d = 4 
with Z = and d = 2 with / = 0, 1, and 2. 



We address d = 4 with / = first. The sum (4.291 contains the single term corresponding 
to k = 0; thus the singular terms are 



-2 /-A 



and the logarithmic terms 



2 A ( 2 A 
ln(az/2) — j (az) - ln(bz/2)Bi (bz) = ln(z) B + 0(z 2 ) ) + 0(1), 

7T \ 7T 



since jo(0) = io(0) = 1. To make i?(z) continuous at r = 0, we must then take 



A=-%B and A = —B. 
2 b 2 2 



Since a < b, this is only possible when both A and B are zero. 

When d = 2, we have s = and ni(z) = Ni(z) and fc;(z) = K\{z). We first consider I = 1 



and 2. In these cases, we only have one singular term, corresponding to fc = in the sum (4.29) 
Assume B ^ 0; as before we may take B = 1. Then the linear combination 



R(z) 



tt /a 
2 



(|) Ni(az) + K t (bz) 



is continuous at z = 0, and so we have by (4.281 that 



R(z)=ln(z) J / (az) + (-l) i ln(z)/ i (^)) +0(2 
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If I = 1, we have 



R'{z) = ln(z) aJ[(az) - Mf(te)) + 0(1) 

= \M*) (y-& + 0(z 2 ))+0(l). 

Since a < b, we have — > oo as z — > 0. Thus we must take B = (and hence A — in 
order for -R(z) to be continuous with R'(z) continuous at z = 0. 
If Z = 2, we have 



and so 



R(z) = ln(z) (j^Y Maz) + I 2 (bzfj + 0(z 2 ), 



*,^ 1 » W ((^ 2 < 1 ^, + « ; (M) + o(,), 



and 



i?"(z) =ln(z) ^) 2 a 2 4'(az)+6 2 ^'(6z)^ +0(1). 

Because J% (0) = 7^' (0) = 1/4, we have that as z — > 0, R"(z) behaves like 

ln(^(Q 2 « 2 + & 2 )+0(l), 

and so we find |i?"(z)| —> oo as z — > 0. Thus we must take A and -B both zero in order for R(z) 
to be continuous with R" (z) continuous. 



When / = 0, the sum (4.29) is in fact empty, and the logarithmic terms equal 



In : ( — Jo(az) + BI (bz) j f ()(] i 



Thus R(z) = AN (az) + BK (bz) is continuous at z = if 



A = - n B. 
2 
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As before, we assume B = 1. Then 



R(z) = - ln(z) ( JqM - I (bz)) + 0(1). 

Thus the first derivative is 

R'(z) = - ln(z) (aJ^az) - bl' {bz)) + 0(z) 

and the second derivative is given by 

R"{z) = -]n(*) (a 2 J^'(az) - b 2 I'^bz)) + 0(1) 

Then as z — > 0, we see |i?"(z)| — > oo. Thus we must take ^4 and i? both zero in order for R(z) 
to be continuous with R"{z) continuous. □ 
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CHAPTER 5 
The unit ball 



We will use the eigenf unctions of the ball as our trial functions in the proof of the isoperimetric 
inequality, Theorem [T] Happily, the full set of solutions for the ball can be found exactly in 
terms of Bessel and modified Bessel functions, and we can in fact identify the fundamental mode. 
In particular, the fundamental mode will be proved to have angular dependence. 

We will focus on the unit ball, since the solution of our eigenvalue problem for any ball can 
then be obtained by scaling. We will show in Theorem [2] that all eigenfunctions will be of the 
form Ri(r)Yi(9), where Ri is a linear combination (depending on r) of ultraspherical Bessel and 
modified Bessel functions of order I, and Y\ is a spherical harmonic. 

Spherical harmonics 

In the case where SI is the ball, it is natural to consider spherical coordinates. Let r be the 
radius and 9 be the remaining angular information. Consider Laplace's equation A/ = 0, with 
/ a function on K d . The Laplacian can be written in spherical coordinates as 

a ^ i-L 1. 

A = d„ H —o r + ^Ag, 

where we give the name As to the angular part of the Laplacian. Separating variables so that 
/ = R(r)Y(§), we obtain 

R" + ^—±R> - l ^ l + d ~ 2) R = and A S Y = -1(1 + d - 2)Y. 

Using our earlier notation for the surface Laplacian, we have Aqq = ^Ag, when f2 is the ball 
of radius R. The parameter I appearing in the separation constant 1(1 + d — 2) must be an 
nonnegative integer in order for solutions to exist. The solutions to A$Y = —1(1 + d — 2)Y are 
called the spherical harmonics. For each I, we choose a spanning set {Y{\ of such solutions that 
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are orthonormal with respect to the L 2 (dM d ) norm. Because the eigenvalues are real, the Y\ 
may be chosen to be real-valued. However, they are traditionally chosen to be complex-valued, 
and so will be treated as possibly such in the proof of this chapter's main result, Theorem [2j 



Factoring the eigenfunction equation 

Proposition 8. Let r > and uj be any positive eigenvalue of the free plate when fl is the unit 
ball. Then the corresponding eigenfunctions in Theorem^can be written in the form R{r)Yi{9), 
where Yj is a spherical harmonic of some integer order I and 

R{r) = ji(ar) + ~fii{br), 

where a and b are positive constants depending on r and ui as follows: b 2 — a 2 = r and a 2 b 2 = uj, 
and 7 is a real constant given by 

-a 2 jf(a) 
7 b 2 i'{{b) ' 

Proof. We first show that eigenfunctions can be written as a product of a radial function with 
a spherical harmonic, and then give the exact form of the radial part. 

Write A := A 2 — tA. By Proposition [l] each eigenvalue ui has finite multiplicity, and so 
the corresponding space of eigenfunctions X u is finite-dimensional. Because As is independent 
of r, it commutes with the Laplacian A and hence with our operator A :— A 2 — rA. Thus 
Ag maps into itself. The operator Ag is symmetric, and so diagonalizable on the finite- 
dimensional space X u . The eigenfunctions of Ag on dW d are the spherical harmonics; on M d the 
eigenfunctions have the form R{r)Yi{9). Thus we can choose our eigenfunctions of A to have 
this form. That is, A and As are simultaneously diagonalizable. 



To find the precise form of R, we factor the eigenvalue equation (1.2), obtaining 



(A + a 2 ){A - b 2 )u = 0, (5.1) 

where a and b are positive real numbers satisfying b 2 = a 2 + t and lo = a 2 (a 2 + r). That is, 
a 2 = W (t/2) 2 + uj — t/2 and b 2 = W (t/2) 2 + u> + r/2. The eigenfunctions u will then be linear 
combinations of the solutions v and w of each factor: 

(A + a 2 )v = and (A - b 2 )w = 0. (5.2) 
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Each of these is separable in spherical coordinates, with angular equation 

A$Y = —l(l + d— 2)Y for some nonnegative integer I. The radial equation for v is a rescaling of 



the ultraspherical Bessel equation (4.1 ) with order I and the radial equation for ui is a rescaling 



of the ultraspherical modified Bessel equation (4.2) with order I, hence 



v = (Aj h (ar) + Bn h (ar)^Y h and w = (ci h (br) + Dk h (&r)) Y h , 

for some nonnegative integers li, I2 and real constants A, B, C, and D. From the diagonalization 
argument above, we know u = i?(r)Yj(0), so all the orders must agree: I = l\ = ?2- Thus 



solutions of the eigenvalue equation (1.2) have the form 



u(x) - R{r)T{9) = [Aji{ar) + Bn t (ar) + C H {br) + Dk h (br) }Yi{6) . 

However, we have from Proposition [2] that the eigenfunctions are smooth on M d . The spherical 
harmonics Yi have no radial dependence; thus we must have the radial part R(r) be smooth 
for r € [0,oo). When I = 0, the spherical harmonic Yq is constant, and we must also require 
R'(0) = in order for u to be smooth. When 1 = 1, the spherical harmonics Y\ can be given by 
Xi/r, where Xj are the coordinate functions. Then along the x^-axis, R{r)Yi{9) = R(r)xi/r = 
R(r) sign(:ri). This function is continuous at the origin only if R(r) vanishes at r = 0. By 
Lemma |4.7[ there is no nontrivial linear combinination of Bessel functions of the second kind 
which satisfies these conditions; thus B and D are both zero. Denote CJA by the constant 7; 
then we have 

u{x) = R{r)T{6) = (ji(or)+7»i(&r)W0)- 



The constant 7 must be chosen so that u satisfies the natural boundary condition u rr = at 
1 (see (2.4) in Proposition [6| ; hence we have 



r 



7 bH'{{b) ' 

and so 7 is real- valued. □ 

The fundamental mode of the ball 

In this section, we identify the fundamental mode of the ball for positive tension, proving: 
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Theorem 2. For r > 0, the fundamental mode of the unit ball has angular dependence, and 
has the form 

ui(r,0) = (ji(ar)+7«i(&r))*i(0), 

with a, b, 7 real constants, with a and b positive and depending on t and uj as follows: b 2 — a 2 = t 
and a 2 b 2 = uj, and 7 given by 

,, = -° 2 J?(a) 
7 t?i»{b) ■ 

Thus in dimension 2, 

sin(0) 



til {r,6)=(j 1 (ar)+ 1 I 1 {br)) I 



cos(0) 



Proof. The proof will have two parts. First we show that for any radial function R(r), the 
Rayleigh quotient Q[RYr] is minimized when 1 = 1, among all / > 1. Then we show that of all 
nonconstant eigenstates with I = and 1 = 1, the lowest eigenvalue corresponds to I = 1. Note 
that when I = 0, the spherical harmonic Y is the constant function, and so I — corresponds 
to purely radial modes. 

[Part 1.] We will show that for any fixed smooth radial function R, the Rayleigh quotient 
Q[RYi] is an increasing function in I for all I > 1. Then by the variational characterization of 
eigenvalues, we see that the lowest eigenvalue corresponding to an eigenfunction with angular 
dependence (i.e., I > 1) occurs when 1 = 1. 

Considering the numerator and denominator separately, we will use the L 2 -orthonormality 
of the spherical harmonics to simplify the angular parts of the integrals. 

The denominator of our Rayleigh quotient is, for u — RYi, 



f iRY^dx^ f . 
Js d Jo 



R 2 r 4 ' 1 dr, 



and so is independent of I. So it suffices to show that the numerator is an increasing function 
of I for I > 1. 

Recall the numerator of the Rayleigh Quotient is 



N[u] = [ \D 2 u\ 2 + r\Du\ 2 dx. 
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We use the pointwise identity of Fact [2] to rewrite the Hessian term as 



\D 2 u\ 2 = * (A(|£>u| 2 ) - D(Au) ■ Du - D(Au) ■ Di 



(5.3) 



Because our region fl is the unit ball, we may use spherical coordinates, noting §^ — u r . Recall 
that 



Du = u r r H — Vsit, 
r 



(5.4) 



where f = x/r is the unit normal. Recall from the derivation of the boundary conditions on the 
ball from Chapter [2] that ^Vs = grad af2 is the surface gradient and ^Ag = divgngrad ao is the 
Laplacian on the boundary of the ball. 

Note by the Divergence Theorem on dM d , we have for any function /, 



/ \&sfdS=[ div aB d ( grad B£i f) dS = 0. 

JdB d r JdB d K ' 



(5.5) 



Further exploiting orthonormality of the Y[, we see that 



f \V s Y l \ 2 dS = l{l + d-2)=:k. 



(5.6) 



Thus when u = RYi, we can rewrite (5.3) as follows: 



- / (A{\Du\ 2 ) - D(Au) ■ Du - D{Au) ■ Du) dx 
2 Ja d ^ ' 



^_ / d-1 k \ „ /_ d-l_ k_. . 
— Du ■ D u rr H u r ~u ~ Du ■ D u rr H u r ^u ) dx 



(noting A s Yi = -kY t ) 

1 r ( d 2 d-1 d 

dr 2 r dr 



[tifW + ^sYA 2 



— Du ■ D I u r 



d-l 



-u r 2 U ) — Du ■ D \ u r 



d-l_ k , 

U r 2 U I 



with this last by noting that (5.5) gives us 



L ^ As ( (i? ' )2|y * 

JR d r \ 
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Expanding the integrands, we see (5.3) becomes 



Mr (R 



?n2 

"tRR' 



\WsY^[^-XRR' + k - d + A R^dS 



Then integrating the above over M d using (5.6) and the orthonormality of the Yi, we obtain 



\D 2 u\ 2 dx 



o 



(i?" 



2k + d-l.,. 2 6fc . , k(k-d + 4) u2 \ , 1 



'-(my - —~rr 



R z r d - L dr 



(R"f + ^(R') 2 + ^(rR'-lR 



k(k-d-l/2) R2 \ d _ ldr 



with this last equality by completing the square. 

We now examine the gradient term in the numerator of the Rayleigh quotient: 



\Du\ 2 dx = 



1 



u r \- + ^|V s u| ) dx 



i| 2 + ^W*| 2 ) dx 



[ ((R') 2 \Yi 



(R 



./\2 



R 2 r d ~ l dr, 



again by (5.6) 



Combining these results, the numerator of the Rayleigh quotient can be now written with 
all /c-dependence (and hence Z-dependence) explicit: 

f 1 (2k f , 3 \ 2 k(k-d- 1/2) 2 kR 2 \ d , , 
N[u] = J I— [rm--R] + -i j '—tR 2 +T— \ r" l ,lr 



- I | (R") 2 + d -^(my + T(my) ,-■> '<//■ 



Recall r > 0; then the above is increasing with k for k > d + 1/2. Recall fc = Z(Z + d — 2) is 
increasing as a function of Z; then all terms involving k are increasing functions of Z for I > 2. If 
Z = 1, the expression k(k — d — 1/2) becomes — 3(cZ — l)/2, which is negative for all dimensions 
under consideration. If I = 2, we find k(k — d — 1/2) = 2d(d — 1/2) > 0. Thus each term 
involving Z is increasing as a function of Z for all Z > 1 . 

Thus for any fixed radial function R, the numerator _ZV[i?Y;] is an increasing function of Z for 
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I > 1. 

[Part 2. J We now show that the lowest eigenvalue corresponding to an eigenfunction of the 
form ui = (ji(ar) + "fii(br)^Yi with I = 1 is less than the lowest positive eigenvalue for I = 0. 
Recall b — \] o? + r. Let aoo denote the first positive zero of j'i(a). All eigenfunctions u 



satisfy the natural boundary conditions, Mu = and Vu = on dfl, as given in (2.4) and (2.5 1 
for the ball. Since all eigenfunctions are linear combinations of ji(ar)Yi(0) and ii(br)Yi(9), we 
must have some nontrivial linear combination satisfy the homogeneous linear equations 




Thus we need the determinant 



Wi(a) := dct 



Mji(ar) Mii(br) 
Vj^ar) Vk(br) 



to vanish. Using the form of natural boundary conditions for the ball given in Proposition [6j 
we have 

Mji(ar) = a 2 j"(ar) and Mi t (ar) = b 2 i'{{br). (5.7) 
The ji(ar) and ii(br) are rescaled ultraspherical Bessel and modified Bessel functions, so by the 



factorization (5.2), we have 



Aji(ar)Yi{6) = -a 2 ji(ar)Yi(e) and Ai,(6r)l!(0) = b\{br)Yi{6). 



Then noting r = 1 on dM d , the "V" boundary condition terms from in (2.5) can be rewritten 
as follows: 



Vji(a) = raj; (a) + k(aj[(a) - j/(a)J + a 3 j z '(a) 
Vii(b) - rbi\{b) + k(bi[{b) - i,(6)) - 6 3 zJ(6) 
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Combining the above with (5.7) and substituting r = b 2 — a 2 , we find 



Wi{a) = a 2 jf (<*)( ~ a 2 U\{b) + k(bi[{b) ~ 

-b^'/^f ab^Ka) + k(af l (a)- Ma)) 



(5.8) 



Given r, the roots of Wi(a) will determine the eigenvalues by the relation u = a 2 (a 2 + t). 
The parameter r is positive, so w increases with a. Note that if a — is a root, it corresponds 
to u> = 0. 

Therefore to show that the lowest nonzero eigenvalue corresponds to I = 1 and not ^ = 0, we 
show that the first nonzero root of W\(a) is less than the first nonzero root of Wo (a). 
First we consider I — 0. Here k — 0, so we look for solutions to: 



= W a (a) 

= a 2 Jo »( - aX(fe)) - bX(b)(ab 2 f (a) 
= a 4 bj' 1 {a)i 1 (b) + abH' 1 (b)j 1 (a), 



by (4.4 1 and (4.7). The functions io(b) and Iq(6) are positive for b > 0, as noted in Chapter [4] 



by the power series expansion (4.12). Similarly, ji(a) and j[(a) are positive on (0,aoo) by 



Lemma 4.2 so Wq(ci) > on (0, a c 



Now consider 1 = 1. The constant k = d — 1, so we have 



Wi(o) = a 2 ^'(«)( - « 2 KW + (d- l)(Wi(&) - H{b)) 
- b 2 i'((b)(ab 2 j[(a) + (d- l)(aiJ(o) - Ji(o)) 
= a 2 6ji'(a)( - a 2 «i(6) + (d- l)i 2 (&)) 



a6 2 



ii'(6)(6 2 il(a)-(rf-l)i 2 (a)), 



(5.9) 



(5.10) 



by the Bessel identities (4.4) and (4.7). As a — > 0, the first term in (5.10) behaves like 



a 2 V^f{(a)((d-lMV^)) 



and so its sign is determined by that of j"(a). By Lemma 4.4 j"(a) is negative for all a G (0, a c 
Hence the first term of Wi(a) is negative near a — 0. 
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Similarly, as a — » 0, we see .72(a) —> 0, and so the second term behaves like 



which is negative by Lemma 4.2 Therefore Wx(a) < near a = 0. 

At a = a,oo, we have Ji(aoo) = by definition of doo! also note Ji(aoo) > and j"(a,oo) < 

r. By Proposition [TJ we have > d > d — 1 for 



by Lemma 



4.1 



and 



4.4 



Write b Q 



d > 3; for d = 2, we have « 1.84 2 > d - 1. Thus a^, - (d - 1) > 0, and so (5.9 1 gives 



Wi{a O0 ) = -aLj'i(aoo)(feoo(aL - (d- l))t'i(&oo) + (d - l)*i (&«>)) 



+ &ooh ( & oo)(d - l)ji(aoo)- 



Hence both terms in Wi(aoo) and since Wi is continuous, it must have a zero in (0, doo)- Thus 
the lowest nonzero root of Wi(a) occurs when 1 = 1, not I = 0, and so the lowest eigenvalue 
uj = a(a 2 + t) occurs when 1 = 1, □ 

We have now assembled all of the tools we will need to tackle the proof of the main theorem. 
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CHAPTER 6 

Proof of the free plate isoperimetric inequality 



In this chapter, we prove Theorem [T] 

Among all regions of a fixed volume, when r > the fundamental tone of the free plate is 
maximal for Q a ball. That is, 

wi(fi) < Wi(f2*), with equality if and only if is a ball. (6-1) 

For simplicity, in this section we will write uj instead of lu\ for the fundamental tone of the 
free plate with shape fl; the fundamental tone of the unit ball will be denoted by u)*. When 
dependence on the region and the tension r need be made explicit, we write w(r, f2) for 
the fundamental tone and Q T ,n for the Rayleigh quotient. We shall also need the notation 
sf2 := {x £ M. d : x/s £ £1} for s > 0. Note also that r > throughout this section. 

The proof of Theorem [T] will proceed from a series of lemmas, following roughly this outline: 

• Scaling 

• Definition of trial functions 

• Concavity of the radial part of the trial function 

• Evaluation of the Rayleigh Quotient 

• Partial monotonicity of the integrand in the numerator and denominator 

• A collection of lemmas needed to establish this monotonicity 

• Rearrangement 

• Proof of the theorem. 

We begin by rescaling to reduce to domains having the same volume as the unit ball 



l d (Lemma 6.1). Adapting Weinberger's approach for the membrane [52], we construct in 
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Lemma 6.2 trial functions with radial part p matching the radial part of the fundamental mode 



of the ball. We follow by proving in Lemma 6.3 a concavity property of p that will be needed 



later on. We next bound the eigenvalue w by a quotient of integrals over our region O, both 



of whose integrands are radial functions (Lemma 6.4 1. These integrands will be shown to have 
a "partial monotonicity" . The denominator's integrand is increasing by Lemma |6.5| and the 
numerator's integrand satisfies a decreasing partial monotonicity condition by Lemma |6.6| The 
proof of Lemma |6.6| becomes rather involved and so is broken into two cases, Lemma |6.7| for 
large r values, and Lemma |6.8| for small values of r. The latter in turn requires some facts about 
particular polynomials, proved in Lemmas |6.9| and |6.10| We then exploit partial monotonicity 
to see that the quotient of integrals is bounded above by the quotient of the same integrals 
taken over fl* , by Lemma |6.11| Finally, we conclude that the quotient of integrals on ft* is in 
fact equal to the eigenvalue uj* of the unit ball. From there we deduce the theorem. 
Our first lemma is a scaling argument. 

Lemma 6.1. (Scaling) For all s > 0, we have 

oj(t, Q) — s 4 oj(s~ 2 t, sfl). 

Proof. For any u £ H 2 (Jl) with J^udx — 0, let u(x) = u(x/s). Then u is a valid trial function 
on sf2 and so 

J sn \D 2 u\ 2 + s- 2 r\Du\ 2 dx 
_ Ln l s " 2 (- c,2 '")0 c / s )| 2 + fl _a r|«- 1 (I>«)(x/a)| a da! 



s- 4+d J n \D 2 u\ 2 + T \Du\ 2 dy 



J sn u{x/s) 2 dx 
r)n\ 2 dn 

taking y = x/s, 



s d J„ u 2 dy 

Now the lemma follows from the variational characterization of the fundamental tone. □ 



Once we have established inequality (6.1) for all regions f2 of volume equal to that of the 



unit ball and all r > 0, we obtain (6.1 1 for regions of arbitrary volume, since 



lj(t,Q) = s 4 uj{s- 2 T 1 sn) < s 4 uj(s- 2 T,sn*) =cj(t,Q*), 
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for all s > 0. 

Next, inspired by Weinberger's proof for the membrane |52j . we choose appropriate trial 
functions from the fundamental modes of the unit ball. In the following lemmas, we take 



R{r) = ji(ar) + jh(br) 

to be the radial part of the fundamental mode of the unit ball. Recall a and b are positive con- 
stants determined by r and the boundary conditions, as in the proof of Theorem 2 in Chapter [5j 
The constant 7 is positive and determined by a, t, and the boundary conditions to be 

7«^>0. (6.2) 
Recall also that R(r) > on (0, 1] and R'(l) > 0. 

Lemma 6.2. (Trial functions) Let the radial function p be given by the function R, extended 
linearly. That is, 

R(r) when < r < 1, 

R(l) + (r - l)-R'(l) whenr>l. 

After translating Q suitably, the functions Uk — Xkp(r)/r, for k = l,...,d, are valid trial 
functions for the fundamental tone. 

Proof. To be valid trial functions, the Uk must be in H 2 (Q). Because Vl is bounded in M. d , 
the only possible issue is at the origin. The series expansions ( |4.11 1 and ( 4.12[ ) of j\ and i\ 



p(r) = 



respectively give us that R(r)/r approaches a constant as r — > 0. thus, Uk £ H 2 (£l) as desired. 
They must also be perpendicular to the constant, that is, 

^ - k dx = for k = 1, . . . ,d. 

Q r 

We use the Brouwer Fixed Point Theorem to translate our region so that the above conditions 
are guaranteed; here again we follow Weinberger [52J. Write x = (x\, . . . ,Xd) and consider the 
vector field 

x{v)= f 4^i (x „ u)da , 

Jn F ~ v\ 

The vector field X is continuous by construction. Along the boundary of the convex hull of Q, 
X(v) is inward-pointing, because p > and the entire region lies in a half-space to one side 
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of v. Thus by the Brouwer Fixed Point Theorem, our vector field X vanishes at some v in the 
convex hull of f2. If we first translate f2 by v, then we have X(0) — f n p(r)x/r dx = 0. This 
gives us f n Uk dx = 0, as desired. □ 

We will need one further fact about our radial function p. 

Lemma 6.3. (Concavity) The function p"(r) < for r g [0,1], with equality only at the 
endpoints. 

Proof. First note that on [0, 1], the function p = R. We see 

R"{r)=a 2 j'{{ar)+ 1 bH'{{br), 



which is zero at r — because the individual Bessel derivatives vanish there, by the series 



expansions in the proof of Lemma 4.6 At r — 1, the function R" vanishes because of the 
boundary condition Mu = 0. 

The fourth derivative of R is given by 

R""{r)=a i j1"{ar)+ 1 bH'{"{br). 



Because all derivatives of ii(z) are positive when z > 0, the second term above is positive 



on (0, oo). We have by Lemma 4.5 that j""(z) is positive on (0, aoo]. Thus R""{r) > on (0, 1], 
and so R"(r) is a strictly convex function on [0,1]. Since R" = at r = and r = 1, the 
function R" must be negative on the interior of the interval [0, 1]. □ 

We now bound our fundamental tone above by a quotient of integrals whose integrands are 
radial functions. Write 

, r r 1 / r/\2 i 3(a— 1) /2 , 2 T(d — 1) 2 
N \p] ■= (P ) + i — KP-rp) +T{p ) + 2 — ■ 



Lemma 6.4. (Evaluation) For any Q, translated as in Lemma 6.2 we have 



L N[p] dx , 
w < J " 2^ (6.3) 



toitt equality if Q = £1* 
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Proof. For u k denned as in Lemma 6.2 we have 



, S a \D 2 n t \ 2 + ADu t fix 

" SQW = ' 

from the Rayleigh-Ritz characterization. We have equality when Q — Q* because the u k are the 
eigenfunctions for the ball associated to the fundamental tone by Chapter [5] Multiplying both 
sides by J n \uk\ 2 dx and summing over all k, we obtain 

. d 



f J2\u k \ 2 dx< j Y,\D 2 u k \ 2 +TY,\Du k \ 2 dx (6.4) 
"' Q k-l ^ n k=l k=l 



again with equality if ft — f2*. From Appendix [j Fact [j| we see inequality (6.4) becomes 



-• / /-/.- JJ »" r + - rp') 2 + r{p') 2 + U ,r j ,/.,-. 



once more with equality if Q is the ball 0*. Dividing both sides by p dx, we obtain (6.3). □ 



We want to show the quotient in Lemma |6.4| has a sort of monotonicity with respect to the 
region J7, and so we examine the integrands of the numerator and denominator separately. 

Lemma 6.5. (Monotonicity in the denominator) The function p(r) 2 is strictly increasing. 

Proof. Differentiating, we see 



p'(r) 



j'i(ar) + ji'i(br) when < r < 1, 
R'(l) whenr>l. 



Obviously i'\(br) > 0. Because we have a < a x from Chapter[5j the function j'i(ar) is positive on 
[0, 1]. Thus p'{r) is positive everywhere, and p (and therefore p 2 ) is an increasing function. □ 

The monotonicity result for the numerator is rather more complicated and requires several 
lemmas. We do not need to prove the integrand of the numerator is strictly decreasing; a weaker 
"partial monotonicity" condition is sufficient. We will say a function F is partially monotonic 
for f2 if it satisfies 

F{x) > F(y) for all x £ fi and y £ O. (6.5) 
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Lemma 6.6. (Partial monotonicity in the numerator) The function 



N[P] = {P"Y + 



3(d- 1) 



(p-rpiy + r (piy + id-i) 



satisfies condition (6.5) for f2 the unit ball 



Proof. Given that p" < on (0, 1) and equals zero elsewhere by Lemma 6.3 the function (p") 2 
satisfies condition (6.5) for the unit ball. The derivative of the function r(p') 2 with respect to 



r is 2rp' p" , and hence negative on (0, r) and zero everywhere else. Thus r( j o') 2 is a decreasing 
function of r. It remains to show that the remaining term 



is also a decreasing function of r. Differentiating, we see 



Now, p — rp' — at r = and 



dr 



(P ~rp') = -rp" 



so by Lemma 6.3 (p — rp 1 ) is positive on (0, oo) and vanishes at zero. Thus in order for h(r) to 
be decreasing, we must have 



6 



(p - rp') + 3p" + rp > 0. 



(6.6) 



Let A rj o := p" — (d — l)r (p — rp'). Recall from the Bessel equations (4.1 1 and (4.2 ) that 



A r ji(ar) = —a 2 ji(ar) and A r ii(br) = b 2 i\{br). 



(6.7) 



Then on the interval [0, 1] 



^(P- rp / ) + 3p" + rp= ^(p-rp') + 3f A r p + — —{p-rp') \ +rp 

= - r P') + 3 ( - « 2 JiM + 7^ 2 n(&r)) + 



TP, 



with the last equality by (6.7). Considering the first term of the last line above, we see by (4.4) 
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and (4.7), 



^(p-rp') = ~(arj 2 (ar) - brji 2 (br)^ 

= {a 2 [ii{ar) + j 3 (arf) + ^b 2 (i 3 (br) - i x {br) 



with the second equality by (4.3), (4.6), and simplifying 



Therefore our quantity of interest in (6.6) can be bounded below in terms of j{s and iis: 



(> 



{p -rp')+ 3p" + rp 



3a 2 \ . . , 3a 2 (d+l) . , , 
T - I Ji(ar) H j— - — 33(ar) 



d+2 



+ 7 t + 



3b 2 
d + 2 



h(br) + 7 



d + 2 
3b 2 (d + l) 



d + 2 



i 3 (br 



3a 2 (d+l) . . , 36 2 (d+l). , 
> y —^J 3 (ar)+ 1 ; ■ - J i 3 {br) 



d + 2 



3a 2 
d+2 



+ 7 t + 



d + 2 
3b 2 
d + 2 



with the inequality by ji{ar) < ii{ar) < ii(br), since r > and so a < b. 

The function i 3 is everywhere positive. We have a < a^, so ar < on [0,1] and hence 



the functions j 3 (ar) and ji(ar) are positive on [0, 1] by Lemma 4.1 The remaining factor is 



positive for all r > by Lemmas 6.7 and |6.8| (to follow), thus establishing (6.6) and completing 
the proof. □ 

We establish the positivity of the remaining factor first for those values of 
t > 9/(d + 5); the proof for the remaining r values is more complicated and is treated in another 
lemma. 



Lemma 6.7. (Large t) We have 



for allr> 9/(d + 5). 



3a 2 n 

T - > 

d + 2 



(6.8) 



Proof. We use the bounds we established for w(r) in Chapter [2j 

Recall that the first free membrane eigenvalue for the ball is p\ = a 2 ^. Lemma 3.1 and 
Proposition [7] together give (d + 2)r > ui* > rd. Because lo* — a 4 + a 2 r from our factoring of 
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the eigenvalue equation in Chapter [5j we obtain inequalities relating r and a: 



4 4 

a a 
d-a 2 >T> d+2 -a 2 ' 



(6.9) 



with the upper bound holding only if a 2 < d. 
Using the lower bound, we see 



3a 2 a 

T - ~ > 



3a 2 



rf + 2' d + 2-a 2 d + 2 

_ a 4 (rf + 5) -3a 2 (rf + 2) 
~~ (d + 2){d + 2- a 2 ) 

which is nonnegative whenever a 2 > 3(d + 2)/(rf + 5). When a 2 < 3(d + 2) /(d + 5), we have 



3a 2 9 

T -d-^-2 >T -dTl > ^ 



by our choice of t. 



Lemma 6.8. (Small r) We have 



□ 



3a 2 / 36 2 . 

t ~ -. r + 7 r + -^-^ ) > 



d + 2 



d + 2 



(6.10) 



for all0<r< 9/(d + 5). 

Proof. The proof will proceed as follows. For < t < 9/ (d+5), we restate the desired inequality 



(6.10) as a condition on 7, (6.12). We then use properties of Bessel functions to establish a lower 



bound on 7 in terms of a rational function of a; we then show this function satisfies (6.12). We 



will need to treat the cases of d > 3 and d = 2 separately, because the two-dimensional case 
requires better bounds than we can derive for general d. 



First note that b 2 = a 2 + t, so the inequality (6.10) is equivalent to 



T > 



3a 2 (l- 7 ) 



(d + 5)7+ {d + 2)' 



(6.11) 



Using the lower bound on r in (6.9), we see that the above will hold if 



3(d + 2)-a 2 (d + 5) , 
7 " (3 + a 2 )(d + 2) ' 7 



(6.12) 
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We need only show that (6.12) holds for all < t < 9/(d+5). We will use Taylor polynomial 
estimates to bound 7 below by a rational function. From Lemma |4.6[ we have 



j"(z) < -diz + d 2 z 3 
6 

i'{(z) < d x z + -d 2 z 3 





on [0, v / 3{d + 2)/{d + 5)} 
on [0, V3], 



These bounds apply to z = ar and z = br respectively, when r € [0, 1], as we show below by 
obtaining bounds on a 2 and b 2 . 



To derive our bound on a 2 , we note that the lower bound of (6.9) together with our assump- 
tion t < 9/(d + 5) implies 



< 



so that 



d + 2-a 2 d + 5' 



(d + 5)a 4 + 9a 2 - 9(d + 2) < 0. 



The lefthand side is increasing with respect to a 2 and equals zero when a 2 = 3(d + 2)/(d + 5). 
Hence a 2 < 3(d + 2)/{d + 5) and the bound on j"(z) holds for z = ar when t < 9/(d + 5). We 
use these to obtain a further bound: 



> T 



> 



9 



d + 5 

„4 



d + 2- a 2 d + 5 

a 2 + 3 / 2 3(rf + 2) 
d + 2 -a 2 \ d + 5 



and so we have a 2 < d. 

To bound b 2 , we use b 2 = a 2 + r and obtain 



l2 2 ^ 3 d + 2 9 
b 2 = a 2 + t < — -!- + - — 
~ d+5 d+ 



= 3, 



and so b 2 < 3. 



We also have, from (6.9), 



da 2 , 9 (d + 2)a 2 

> b > — 

d-a? d + 2 -a 2 ' 



(6.13) 



with the upper bound holding in this regime because a 2 < d. 
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We also need the following binomial estimate: 



l--x<(l-xf' 2 
2 v ; 



for < x < 1. 



(6.14) 



Using these bounds, we see 



a 2 jf» 



> 



> 



b 2 i'{(b) 

a 2 {d\a — d 2 a 3 ) 
b 2 (d 1 b+(6/5)d 2 b 3 ) 

a 3 (di — d 2 a 2 ) 



( 



3/2 
7 2\ 3/2 



^) (d 1 + (6/5)d 2 ^) 

(d-a 2 )(l- Cl a 2 ) 



> 1- 



d J (d-a 2 + {6/5) Cl da 2 ) 
3a 2 \ (d-a 2 )(6(d + 4)-5a 2 ) 



2d 



(6d(d + 4) - 24a 2 ) 



noting that a 2 /d < 1 and a 2 < 3(d + 2)/(d + 5). 
Thus we have 7 — 7* > if 



by definition (6.2 ) 



by Lemma 4.6 



by (6.131 



writing c\ — d^/di — 5/6(d + 4) 



by (6.141, 



! _ 3o^\ (d-a 2 )(6(d + 4)-5a 2 ) _ 3(d + 2) - q 2 (d + 5) 
~^2d~J 6d(d + 4) - 24a 2 (3 + a 2 )(d + 2) " ' 



or, clearing the denominators and writing x = a 2 , if 

(2d - 3x)(d - x) (6(d + 4) - 5x) (3 + x)(d + 2) - 2d(6d(d + 4) - 24a;) (3(d + 2) - cc(d + 5)) > 0. 

The above polynomial is fourth degree in each of d and x and has the root x = 0; because we 
are only interested in its behavior for x € (0, 3(d + 2)/(d + 5)), we may divide by x and work to 
show the resulting polynomial 

P(x, d) = 24d 4 + 60d 3 - 120d 2 - 432d - 40d 3 x - 119d 2 a; - 6da; + 432a; 
+ 43d 2 a; 2 + 113dx 2 + 54a; 2 - 15da; 3 - 30a; 3 



is nonnegative for x <E (0, 3(d + 2) /(d + 5). This claim is addressed in Lemma 6.9 for d > 3. 

For d = 2, the function P(x, 2) is negative on most of our interval of interest [0, 12/7], and 
so we must improve our lower bound on 7. The derivation follows that of inequality (16.91) in the 
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proof of Lemma 6.7 as follows 



By Lemma 3.1 u>* > a^r, where a^ w 1.84118 is the first zero of J x {z). By Chapter [Hj wo 



a t, giving us 



r < 



Using 6 2 = a 2 + r, we obtain also a bound on b 2 : 



b 2 < 



alo-a 2 ' 



Proceeding as before, we deduce 



3a^\ (al-a^M-Ba 2 ) 
1 2a 2 J 36a 2 , + (6a 2 , - 36)a 2 



with the last again from (6.14). So 7 — 7* > if 



I _ (a 2 -a 2 )(36-5a 2 ) _ 12 - 7a 2 > Q 



2a 2 QO J 36a 2 , + (6a 2 , - 36)a 2 12 + 4a 2 



or, setting x — a 2 , if the fourth degree polynomial 



Q(x) = (l- ~) (a 2 , - x)(36 ~ 5x)(12 + Ax) - (36a 2 , + (6a 2 , - 3%) (12 - 7x) 



is positive on [0, 12/7]. This positivity follows from Lemma 6.10 completing our proof. □ 



Lemma 6.9. The polynomial 



P(x, d) = 24d 4 + 60d 3 - 120d 2 - 432d - 40d 3 x - l\M 2 x - 6dx + 432x 
+ 43d 2 x 2 + 113dx 2 + 54x 2 - 15da; 3 - 30x 3 



is nonnegative for all x £ (0, 3(d + 2)/(d + 5)) and integers d > 3. 

Proof. First note that 3(d + 2)/(d + 5) < 3. We bound P below on the interval x € [0,3] 
by taking x = 3 in terms with negative coefficients and taking x = in terms with positive 
coefficients, obtaining 

P(x, d) > 24d 4 - 60d 3 - 477d 2 - 855d - 810 =: g(d). 
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The highest order term is 24c? 4 , and so g is ultimately positive and increasing in d. Note also 
that 



g'(d) = 96d 3 - 180d 2 - 954d - 855 



g"(d) = 288d 2 - 360d - 954. 



The function g"(d) is a quadratic polynomial with positive leading coefficient and roots at 
d w —1.30 and 2.55; thus g'(d) is increasing for all d > 3. We see that <?'(5) = 1875, so g is 
increasing for all d > 5. Finally, <?(7) = 6876, so for all d > 7 we have g(d) > and hence 
P(x, d)>0 for all d > 7 and x € [0, 3]. 

For d = 3, 4, 5, 6, we look at the polynomials Pd(x) = P(x, d) directly to show that Pd(x) > 
on [0, 3(d + 2)/(d + 5)]. Each Pd is a cubic polynomial in x; its first derivative P' d {x) is quadratic 
and so the critical points of Pd(x) can all be found exactly. 

For d = 4, 5, and 6, direct calculations show P' d < on [0, 3(d + 2)/(d+ 5)] and Pd(3(d + 
2)/(d + 5)) > 0, so P d (a:) > on [0, 3(d + 2)/(d + 5)]. 

For d = 3, our interval of interest is [0, 15/8]. We have a critical point 
c w 1.4 G [0, 15/8], with P^ < on [0,c] and P^ > on [c, 15/8]. The critical value P 3 (c) « 79 
is positive, so Ps(x) > on the desired interval [0, 15/8]. □ 

Lemma 6.10. The polynomial 



is positive on [0, 12/7]. 

Proof. As in previous cases, x = is a root of this polynomial, so we examine g(x) := Q{x)/x. 
The derivative t/(x) is a quadratic polynomial, so its roots can be found exactly. We see that g 
has a critical point c sa 1.4 in [0, 12/7], with g' < on [0, c] and g' > on [c, 12/7]. The critical 



Our final lemma is a simple observation about integrals of monotone and partially monotone 
functions, which is a special case of more general rearrangement inequalities (see [291 Chapter 




value g(c) w 177.8 is positive, so g > on [0, 12/7 



□ 



3]). 
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Lemma 6.11. For any radial function function F(r) that satisfies the partial monotonicity 



condition (6.5) for Cl 



Fdx < Fdx 
n io* 



with equality if and only if Cl — Cl* . For any strictly increasing radial function F(r), 



Fdx> Fdx 
n Jn* 



with equality if and only Cl = Cl* 



Proof. Note that \Q\ = \Cl*\ with \Q\Q*\ = |0*\fl|. Suppose F satisfies for Cl*. The result 
follows from decomposing the domain: 



Fdx= Fdx+ Fdx 

<> Jnnn* Jn\n* 



< [ Fdx+ sup |F(a;)||fi\n*| 

< [ Fdx+ inf |F(af)||n*\n| 

< Fdx+ Fdx 
Jnnn* Jn*\n 

Fdx. 



since F satisfies (6.5) 



Note that if |f2\fT | > 0, either the second inequality or the third is strict by the strict inequality 



in ( 6.5 ) . If F is strictly increasing, then apply the first part of the Lemma to the function —F. □ 



We can now prove our main result. 



Proof of Theorem^ By rescaling as described after Lemma 6.1 it suffices to prove the theorem 
for fl with volume equal to that of the unit ball, so that CI* is the unit ball. We may also 



translate CI as in Lemma 6.2 which of course leaves the fundamental tone unchanged. Then, 



In N \P\ dx 
In P 2 dx 
< /„. N[p] dx 
I n ,p 2 dx 



by Lemma [6^4 



by Lemmas 6.5 6.6 and 6.11 



by applying the equality condition in Lemma |6.4| Finally, if equality holds, then CI must be a 



58 



ball, by the equality statement in Lemma 6.11 



□ 
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CHAPTER 7 

One dimension: the free rod 



The free rod is the one-dimensional case of the free plate. We do not have an isoperimetric 
inequality for the free rod, because all connected domains of the same area are now intervals of 
the same length, identical up to translation. However, we do have a one-dimensional analog of 
Theorem[2j which identifies the fundamental tone of the circular free plate. In higher dimensions, 
we showed the fundamental mode of the free plate under tension had angular dependence; in one 
dimension, we will show the fundamental mode of the free rod under tension is an odd function 
about the center point. In this chapter we prove the one-dimensional analog of Theorem [2] and 
classify most eigenf unctions of the free rod under tension and compression. 

We take f2 = [—1,1]; the general case follows from rescaling and translation. The Rayleigh 
quotient in one dimension is 



Q[u] = 



i £ 1 |u"| 2 + r|u'| a da; 



J_i \u\ 2 dx 

From this quotient, we obtain the eigenvalue equation 



: OJU 



(7.1) 



and the natural boundary conditions 



u" = 
tu' - u'" = 



at x = ±1, 
at x = ±1. 



(7.2) 
(7.3) 



As in the case of the ball in R d , we find the general form of the eigenfunctions by factoring the 



eigenfunction equation (7.1 1. The factorization depends on the sign of the eigenvalue uj and in 



some cases on the value of r relative to lj. We will therefore treat positive, zero, and negative 
eigenvalues as separate cases. In each case, we will use the boundary conditions to help further 
identify form of the eigenfunctions. 



GO 



First we show we need only consider even and odd solutions, before embarking on the classi- 
fication. In higher dimensions, we were concerned with angular dependence in our solutions, so 
it is natural in one dimension to consider the symmetries of the solutions about the center point. 
Note that if u(x) is an eigenfunction, then by symmetry so is u(— x), with the same eigenvalue 
uj. The even and odd parts of u can be expressed as 

. . u(x) — u(— x) , . . u(x)+u(—x) 
u (x) = and u e {x) = , 



which are either solutions of ( 7.1 ) with the same eigenvalue, or (in the case that u is purely odd or 
purely even), one of them is zero everywhere. Because u(x) and u{—x) both satisfy the boundary 
conditions, u e and u Q will also satisfy them. Since every eigenfunction is a linear combination 
of its odd and even parts, it suffices to look only for even and odd eigenfunctions. We will 
refer to eigenvalues associated with even and odd eigenfunctions as even and odd eigenvalues, 
respectively. 



Positive eigenvalues 

In the case of positive eigenvalues, the general forms of the eigenfunctions do not depend on the 
sign of t. As in the case of the plate, positive eigenvalues correspond to vibrational frequencies 
of the rod; the corresponding eigenfunctions are then vibrational modes. In higher dimensions 
and taking fi to be the ball, we saw solutions involving Bessel J/ and /; functions; here the 
corresponding solutions are trigonometric and hyperbolic trigonometric functions. 

Fact 7.1. Given r6l, a positive number uj is an eigenvalue if and only if it has an associated 
eigenfunction u with uj and u satisfying exactly one of the following: 

(i) We have u = u (x) = A sin ax + £> sinh fa; with A, B real nonzero constants satisfying 

B a 2 sin a 
A b 2 sinh b 

and a, b positive numbers satisfying a 2 b 2 — uj, b 2 — a 2 — r, and W (a) = a 3 sin a cosh b — 
b 3 cos a sinh 6 = 0. 
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(ii) We have u = u e (x) = C cos ax cosh bx with C, D real nonzero constants satisfying 

D a 2 cos a 
C b 2 cosh b 

and a, b positive numbers satisfying a 2 b 2 — u>, b 2 — a 2 = r, and W e (a) = a 3 cos a sinh & + 
b 3 sin a cosh 6 = 0. 

We will also prove a theorem identifying the fundamental tone of the free rod under tension:: 

Theorem 2". For all r > 0, we have that the even and odd eigenvalues are interlaced. In 
particular, the fundamental mode is an odd function. 

We prove Fact |7. 1| first. 



Fact \7.1\ We begin by factoring the eigenvalue equation to establish the general form of the 
even and odd solutions. We will then use the boundary conditions to obtain the desired values 
for B/A and D/C and show that they are nonzero and well-defined. 



When u) is positive, regardless of the value of r, we can factor the eigenvalue equation (7.1 1 

as 

{d 2 +a 2 )(d 2 -b 2 )u = Q 1 

where u> = a 2 b 2 and r = b 2 — a 2 . Since oj > 0, a and 6 must be nonzero as well; we take them to 
be positive. Recall that we need only look for even and odd solutions; all others will be linear 
combinations of these. 

Solutions to (d 2 + a 2 )u = are of the form e ±alx ; we express them as the odd and even 
trigonometric functions sinxa and cos ax. Solutions (d 2 — b 2 )u — are of the form e ±bx , or 



sinhba; and cos bx. Any solution u to (7.1 1 is then a linear combination of these terms; since we 



need only consider odd and even functions, we obtain 

u {x) = A sin ax + B sinh bx and u e (x) — C cos ax + D cosh bx. (7-4) 

Applying the boundary conditions to u and u — e, we obtain: 

a 2 ^4sina = & 2 i?sinh6 and a6 2 ylcosa = a 2 6£?cosha (7.5) 

a 2 C cos a = b 2 D cosh b and — ab 2 C sin a = a 2 bD sinh a (7-6) 
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We will now show that for odd solutions u , both A and B are nonzero, and that for 
nonconstant even solutions u e , both C and D are nonzero. 

If A = but B nonzero, then the odd eigenstate is u Q (x) — B smh(bx) and the boundary 



conditions (7.5) become 

b 2 B sinh(b) = and a 2 bB cosh{a) = 0. 
Since b is positive and B is nonzero, the first condition cannot be satisfied. 



If B = but A nonzero, then u {x) = Asm{a) and the boundary conditions (7.5) become 

— Aa 2 sin(a) = and ab 2 Acos(a) = 0. 

Since a, b, and A are nonzero and sin(a) and cos(a) cannot be simultaneously zero, we cannot 
satisfy both conditions. 



If C = 0, our even eigenstate is u e (x) = Dcosh(bx), and the boundary conditions (7.6) 
become 

b 2 D cosh(b) = and a 2 bD sinh(6) = 0. 

Again, b > 0, so the only way both equations can be satisfied is if D = 0, giving us u e (x) = 0a 
constant eigenfunction. 



If D = 0, we have u e (x) = C cos(ax) and the boundary conditions (7.6) become 

a 2 C cos(a) = and — ab 2 C sin(a) = 0. 

Since a, b are nonzero and sin(a) and cos(a) cannot be simultaneously zero, we can only satisfy 
both conditions if C = 0, giving us a constant eigenfunction u e (x) = 0. 

Note that since b > 0, we have sinh6 positive. To obtain the ratios in Fact |7.1| we solve the 



first boundary condition in each of (7.5) and (7.6) 



It remains to show that a and b satisfy the conditions W (a) = for the odd eigenfunction 
u and W e (a) = for the even eigenfunction u e . 

We first consider the odd solution u (x) = Asinax + B smhbx. As in the higher-dimensional 
case, the two boundary conditions Mu = u" = and Vu — tu' q — u'" = give us a pair 
of homogeneous equations which are linear in A and B. Thus we must have the determinant 
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vanish: 

0= (— a 2 sina)(— a 2 b coshb) — (b 2 sinh b)(ab 2 cos a) 
As a and b are nonzero, this simplifies to 

= a 3 sinacoshfr — b 3 cosasinhfe =: W (a). 

Recall b 2 = r + a 2 . Thus, for a given t, our candidate u (x) is a solution to the boundary value 
problem if a satisfies W (a) = and A/B is taken to satisfy either (and hence both) of the 
boundary conditions. 

Now we consider the even solution u e (x) = C cos ax + Dcoshbx. The constants C and D 
are assumed to both be nonzero by Fact |7.1b nd are subject to the constraints (7.6 1. Again wc 
must have the determinant vanish: 

= — a 2 cos a(=Fa 2 6 sinh b) — b 2 cosh bab 2 (=F sin a) . 

That is, 

= a 3 cos a sinh b + b 3 sin a cosh b —: W e (a) . 

Then our candidate u e (x) is a solution to the boundary value problem if a is a root of W e and 
C/D is taken to satisfy either (and hence both) of the boundary conditions. □ 

The proof of Theorem 2" requires the following lemma: 



Lemma 7.2. Let r > and write b = \J a 2 + t. The nontrivial zeros of the functions 

W e (a) = a 3 cos a sinh b + b 3 sin a cosh b and 
W (a) = a 3 sin a cosh b — b 3 cos a sinh b. 

on [0,oo) are distinct and are interlaced. In particular, W has a zero in (0, 7r/2), and the first 
nontrivial zero of W e £ (7r/2,oo). 

Proof. We consider W e first. Recall r = a 2 + b 2 ; viewing a as a variable independent of r, we 
differentiate and obtain 

db —a 
da b 
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Then 



W' e {a) — 3a 2 cos a sinh 6 — a 3 sin a sinh 6 — a 4 6 _1 cos a cosh & 
— 3a& sin a cosh b + b 3 cos a cosh 6 — afo 2 sin a sinh & 
= —a 3 sin a sinh & — Sab sin a cosh b — ab 2 sin a sinh b 
+ (t 2 + 2a 2 r)b~ 1 cos a cosh 6 + 3a 2 cos a sinh b. 

Recall b and hence sinh& are positive. For a in the intervals (tt/2 + 2fc7r,7r + 2fc7r) with fc 
any nonnegative integer, we have sin a > and cos a < 0, giving us W' {a) < 0. When a € 
(37r/2 + 2fc7r, 2-7T + 2kir), the signs are reversed, and W' (a) > 0. 

Similarly, we find that W e (a) > on [2kir, 7r/2 + 2fc-7r] and W e (a) < on [7r + 2fc7r, 37r/2 + 2fc7r]. 
Therefore, We has exactly one zero in each interval (it/2 + kir, (k + 1)tt) for k a nonnegative 
integer, and is nonzero everywhere else 

Now consider W . Differentiating and simplifying, we obtain 

W' a (a) = a 3 cos a cosh b + 3ab cos a sinh +ab 2 cos a cosh b 
3a 2 sin a cosh 6 + (r 2 + 2a 2 r) sin a sinh b. 



Examining the sign of the trigonometric terms as before, we see W (a) > on (2fc7r, it/2 + 2kw) 
and W' (a) < on (ir + 2k?r, 3ir/2 + 2kir) for k a nonnegative integer. We also have W (a) > 
on [tt/2 + 2kn, ir + 2kn] and W a (a) < on [37r/2 + 2kn, 2ir + 2kn}. Thus, W a has exactly one 
zero on (/c7r,7r/2 + kir) for each nonnegative integer k, and is nonzero everywhere else. □ 

We can now prove Theorem 2". 

Proof of Theorem 21' . We see by examining the Rayleigh quotient that for r > 0, all eigenvalues 
are nonnegative. Furthermore, the eigenvalue uj — corresponds to the constant eigenfunction 
and has multiplicity one; all higher eigenvalues are positive. 

By Fact ITU positive eigenvalues are given by uj = a 2 (a 2 + t), with a a root of W e or Wo- It 
is clear uj is increasing in a for r > 0, so the eigenvalues for odd and even modes are interlaced 



if the zeroes of W e and W are. We have this interlacing by Lemma 7.2 Also by this lemma 



W has the first nontrivial root in (0,7r/2), and so the fundamental mode is odd. □ 
As in the higher-dimensional case, Theorem 2" only considers positive tension. If r > 0, the 
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Rayleigh quotient is nonnegative, and so we have only nonnegative eigenvalues. Furthermore, 
the eigenvalue u = is nondegenerate for positive r. 

Our work for negative tension will show that in this regime the eigenvalues often cross, and 
so the theorem does not hold here. 



Zero eigenvalues 

An eigenvalue of zero corresponds to non-vibrational translation of the rod. Other than the 
constant eigenstate u = c, these eigenstates are only possible when r is nonpositive. When 
tension vanishes, we have an odd eigenstate and u) = is of multiplicity 2. For negative tension, 
the rod is under compression, and the eigenvalue u> = is only degenerate when yjrf is a root 
of since or cosx, and in these cases has multiplicity 2. 

Fact 7.3. For all real values of r, uj — is an eigenvalue and the constant function is an 
associated eigenfunction. The eigenvalue uj — is degenerate if and only if it has an associated 
nonconstant eigenfunction u satisfying one of the following: 

(i) (Zero tension) The parameter r = and we have u — u (x) — Ax where A is a nonzero 
real constant. 

(ii) (Negative tension) The parameter r < and u satisfies one of the following: 



(a) We have u — u (x) = Bsinax where a = \/\t\ and B is a nonzero real constant, and 
sin y/\F\ = 0. 

(b) We have u — u e (x) — D cos ax where a — v/]rf and D is a nonzero real constant, 
and cos \/\t\ = 0. 

For all other values of t , the eigenvalue u — is nondegenerate. 

Proof. Note that the form of the solutions depends on the sign of r; we will treat each as a 
separate case. The positive and zero tension cases are fairly easy to solve; we only need to factor 



( |7.1[ ) when t < 0. 

Positive tension. When r > 0, the numerator of the Rayleigh Quotient Q can be zero only 
when u x = almost everywhere. Functions with derivative of zero are themselves constant, and 



so the only solution to (7.1 ) with uj = and r = is the constant solution u — C. Note that in 



this case we have no odd solutions. 
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Zero tension. Here we have U) and r both zero; the eigenvalue equation becomes 



u"" = 0, 

which has general solution 

u(x) = Ax 3 + Bx 2 + Cx + D. 
Our first boundary condition, u"(±l) = 0, gives us 

= u"(±l) = ±6A + 2B, 

so we must have A and B both zero. The second boundary condition then yields u"'(±l) = 0, 
and gives us no additional information. Thus, u{x) = Cx + D is a solution to the boundary 
problem for any C, D G R. The even solutions are then the constant functions u = D; odd 
solutions are given by u = Cx. 

Negative tension. When r < we have a free rod under compression; the proof proceeds 



similarly to the positive eigenvalue case, Fact 7.1 



As in the classification of positive eigenvalues, we factor the eigenfunction equation to get 
the general form of the solution and then apply the boundary conditions to get constraints on 
the coefficients A, B, C, and D, and on possible values for t. 



Since oj is zero, the eigenvalue equation (7.1) factors as 



d 2 x (dl-T)u = 0, 

so we must consider solutions to u" — and u" — tu = 0. Solutions to the first factor are of the 
form u = Ax + B. The solutions of u" ~ tu — are all of the form A cos{ax) + B sin(aa;), with 



a = \/\t\. The even and odd solutions to (7.1 1 are then of the form 

u (x) = Ax + Bsin(ax) and u e (x) = C + Dcos(ax) 

and must satisfy the boundary conditions 

—Ba 2 sin(a) = and At + B(ar + a 3 ) cos(a) = for u D , 

—Da 2 cos(a) = and — D(ar + a 3 ) sin(a) = for u e . 
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Consider the odd solutions first. Because a = y/\r\, we have ar + a 3 = 0, and so the second 
condition for u reduces to A = 0. We then have u (x) = Bsm(ax); we assume B is nonzero so 
that u Q is not even. The first boundary condition then can only be satisfied if a or sin a vanish. 
The parameter r is nonzero, so we must have sin(a) = 0. 

Therefore odd solutions are all of the form Bsinaa; and exist only when \t\ is a root of 
sin y/z. 

Now consider the even solutions. The second condition for u e simplifies by ar + a 3 = to 
= 0, giving us no information. The first condition requires D = or cos(a) = 0. Hence, 
all even solutions will be of the form u e (x) = C + Dcos(ax) with cos(a) = 0. The boundary 
conditions give us no relationship between the coefficients, so C and D may be any real numbers. 
These solutions will only occur only when |r| is a root of sin-^/z. □ 



Negative eigenvalues 

All negative eigenvalues and their associated eigenfunctions correspond to the buckling of the 
free rod under compression. The Rayleigh quotient can only be negative when r is negative, 
which corresponds to a rod under compression rather than tension. 

Finding the solutions for negative eigenvalues is more complicated than either the positive 



or zero eigenvalue cases. The factorization of (7.1 ), and hence the general form of the solutions, 



depends on a relationship between r and w. There are three regimes, named for the behavior 
of the solutions there, and we will prove our results in each regime separately. Unfortunately, 
we cannot complete our classification in the case of the hyperbolic regime; we provide only a 
necessary condition for a number u) to be an eigenvalue in this regime. 

Fact 7.4. Given r , a number u> > — t 2 /4 is a negative eigenvalue if and only if it has an 
associated eigenfunction u with u, uj satisfy exactly one of the following conditions: 

(i) (trigonometric regime) U) > — r 2 /4 and u satisfies one of the following: 

(a) We have u = u (x) = A sin ax + B sin bx with A, B real nonzero constants satisfying 

A 6 2 sin(o) / A acos(b) , . . , \ 
— 'or — = — j—r when sin(a) = I . 



B a 2 sin(a) \ B 6cos(a) 

and a, b positive numbers satisfying —a 2 b 2 = us, —a 2 — b 2 = t, and W (a) 
a 3 sin(a) cos(6) — b 3 sin(6) cos(a). 
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(b) We have u = u e (x) = D cos ax + D cos bx with C, D real nonzero constants satisfying 



C 6 2 cos(6) / C asin(o) , . . 

— = = or — = — — when cosla = 

D a 2 cos(a) \ D osin(a) v ' 

and a, b positive numbers satisfying a 2 b 2 = lu, b 2 —a 2 = t, and W e (a) — 6 3 cos(6) sin(o)- 
a? cos(a) sin(6). 

(ii) (degenerate regime) lu — — r 2 /4 and 

u = u e (x) = C cos ax + Dx sin ax, 



where a = \J\t\/2 rj 1.1394 and is the positive root of sin 2a — 2a/3 = 0, with C, D real 
nonzero constants satisfying 

C 2cos(a) — asin(a) „ 

— = ^ — , . w « -0.4174. 

L> acos(a) 

ffere r « -2.5967 anrf cj ss -1.6856. 

(hyperbolic regime) ^4 number u) < — r 2 /4 is an eigenvalue only if it has an associated eigen- 
function u satisfying one of the following: 

(a) We have u = u (x) — A cos ax sinh bx + B sin ax cosh bx with A, B real constants with B 
nonzero and a, b positive numbers satisfying 

-(a 2 + b 2 ) 2 = lu and 2(6 2 - a 2 ) = r. 

(b) We have u = u e (x) — C sin ax sinh bx + D cos ax cosh ox with C, D real constants with D 
nonzero and a, b positive numbers satisfying 

-{a 2 + 6 2 ) 2 = u and 2(6 2 - a 2 ) = r. 



In the proof of each case, we begin by factoring the eigenvalue equation (7.1), and then 
apply the boundary conditions to obtain conditions on the coefficients A through D, and in the 
degenerate case, on the values of r for which there is a solution. 

The trigonometric regime: — r 2 /4 < u < 

In the trigonometric regime, the eigenfunctions are linear combinations of trigonometric func- 



tions. Here the characteristic function for (7.1) is a fourth-degree polynomial with only purely 
imaginary roots. 
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Proof for the trigonometric regime. In this case, (7.1) may be factored as 



(d 2 x +a 2 )(d 2 x + b 2 )u = 0. 

with a 2 b 2 = — u), —a 2 — b 2 = r, and a and b taken to be nonnegative. Because u) < 0, neither a 
nor b can be zero. The even and odd solutions will then have the forms 



u (x) = Asin(aa;) + B sin(bx) and u e (x) — C cos(ax) + Dcos(bx) 



(7.7) 



with boundary conditions: 



a 2 Asin(a) - b 2 Bsm{b) = 
a 2 Ccos(a) - b 2 Dcos(b) = 



and - ab 2 Acos(a) - a 2 bBcos(b) = 0, (7.8) 
and ab 2 C sin(a) + a 2 bD sin(6) = 0. (7.9) 



If A = 0, then — b 2 Bsxn{b) = and — a 2 bB cos(b) = 0; the only way to satisfy both conditions 
is if B = 0. 

If B = 0, then — a 2 Asin(a) — and — ab 2 Acos(a) = 0, and so A = 0. 

If C = 0, then -b 2 D cos(6) = and a 2 bD sin(fe) = 0, thus D = 0. 

If I? = 0, then — a 2 Ccos(a) = and ab 2 C sin(a) = 0, requiring C = 0. 

For a given even or odd solution, given a and r (and hence b), we must turn to the boundary 
conditions to determine A/B and C/D. It is these ratios that determine the form of the 
solutions, since multiplication by a nonzero constant does not change the eigenvalue or affect 
whether or not the boundary conditions are satisfied. 

If A and B are nonzero, then 



.4 
B 



b 2 sin(6) 
a 2 sin(a) 



or = — 



acos(&) 
b cos (a) 



when sin(a) = 



If C and D are nonzero, then 



C 6 2 cos(6) 



a 2 cos(a) 



or = — 



a sin(6) 
6sin(a) 



when cos(a) = 



Armed with this information, let us consider possible eigenvalue crossings. Given r, suppose 



there exist both u e and u nontrivial solutions as given in (7.7), both with the same associated 
eigenvalue, lu = —a 2 b 2 . Since a and b are completely determined by r and cj, we know A/B and 
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Figure 7.1: Eigenvalue curves with expected crossing points. Green curves are odd eigenvalue 
branches; gold curves are even eigenvalue branches. The purple and blue curves are possible 
odd and even eigenvalue branches, respectively, in the hyperbolic regime. These curves were 
computed under the assumption that all linear combination coefficients are nonzero. 



C/D from above, and so the solutions are determined up to multiplication by a constant. 



By (7.9), we have cos(a) = if and only if cos(6) = 0; similarly sin(a) = if and only if 



sin(b) = by (7.8 1. Now, cos(a) = cos(6) = implies a = (2k + 1)tt/2 and b = (21 + 1)tt/2 for 



some nonnegative integers k, I. Since b 2 



t, solving for r yields 



~((2l + l) 2 + (2k + l) 2 ). 



These values for a and b satisfy the boundary conditions for both the even and odd solutions, 
so the odd and even eigenvalues will coincide at these values of r. Therefore, we expect the 
eigenvalue curves cross at these r values. 

Similarly, sin(a) = sin(6) = implies a — kn and b = lit for positive integers k, I. Again 
we solve for r, obtaining r = — ir 2 (l 2 + k 2 ). Again, these a and b satisfy both sets of boundary 
conditions, and so the eigenvalues would again coincide at these r values. This is illustrated by 
Figure [7] 

The r values for which u> = is degenerate are also illustrated in Figure [7} These de generacies 
occur where even and odd eigenvalue curves cross the line uj = 0, which is the eigenvalue curve 
for the constant eigenfunction. □ 



The degenerate regime: u = — r 2 /4 

In the degenerate regime, the characteristic function has zero as a double root. We will show 
that there is an eigenvalue with an eigenfunction in the degenerate regime for only one value of 
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r; that is, the eigenvalue branches only cross from the trigonometric regime to the hyperbolic 
regime once, when u> ~ —1.6856 and r w —2.5966. The eigenfunction in this case is even. 



Proof for the degenerate regime. In this case, (7.1) factors as 



(d 2 x + a 2 ) 2 u = Q 

With a 2 = — t/2. Since r < 0, a = ^J\t\/2 is strictly positive. By elementary differential 



equation theory, the linearly independent solutions of the differential equation (7.1 ) are cos(aa;), 
sin(aa;), xcos(ax), and xsm(ax). As noted previously, we need only consider even and odd 
linear combinations of these, namely, 



u {x) = Asin(ax) + Bx cos(ax) and u e (x) = C cos(ax) + Dx sin(ax). 



The boundary conditions then yield the relations 

— a 2 Asm(a) — 2a£?sin(a) — a 2 £?cos(a) = 
-a 3 A cos(a) + a 2 B cos(a) + a?B sin(a) = 
— a 2 Ccos(a) + 2aZ?cos(a) — a 2 Dsin(a) = 
a 3 C sin(a) + a 2 D sin(a) - a 3 D cos(a) = 



We will use these to find the values of a (and thus r and w) and the ratios of the coefficients 
that give us solutions to the boundary problems. 

We begin by proving that if u Q (resp u e ) is a nonconstant solution of the boundary problem, 
both A and B (resp. C and D) must be nonzero. 

For the odd eigenfunction, suppose A — but B ^ 0; we seek a contradiction. The boundary 
conditions give us — 2asin(et) — a 2 cos(a) = and a 2 cos(a) + a 3 sin(a) = 0. Since a > 0, we 
obtain — 2sin(a) = acos(a) and cos(a) = — asin(a); thus 2sin(a) = a 2 sin(a). So either 2 = a 2 
or sin(a) = 0. 

If sin(a) = 0, then — a 2 Z?cos(a) = 0. But each term must be nonzero, so this is impossible, 
and a 2 = 2. We still must have cos(a) = — asin(a), but we quickly see this is impossible if 
a = V2. 

Now suppose B — but A ^ 0; we seek a contradiction. Boundary conditions and a > yield 
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sin(a) = = acos(a). Since a^O, both sin(o) and cos(a) must be zero, which is impossible. 

We consider the even eigenfunction next. Suppose D ^ but C = 0; we seek a contradiction. 
The boundary conditions give us 2cos(a) = asin(a) and sin(a) = acos(a), and combining these 
we obtain 2cos(a) = a 2 cos(a). So either 2 = a 2 or cos(a) = 0. 

If cos(a) = 0, then the first boundary condition becomes —a 2 B sin(a) = 0. But a ^ 0, B ^ 0, 
and sin(a) and cos(a) cannot both be zero, so it is impossible to satisfy the boundary conditions 
when cos (a) = 0. 

We must then have a 2 = 2. The boundary conditions still require sin(a) = acos(a), but a 
simple calculation shows that a — y/2 does not satisfy this equality. 

Now suppose C^fl but D = 0; we again seek a contradiction. The boundary conditions 
then yield —a 2 cos(a) = = a 3 sin(a), so then — cos(a) = and asin(a) = 0. Since a is strictly 
positive, we must have both cos(a) = and sin(o) = 0, which is impossible. 

Now that we have shown that A/B and C/D must be defined and nonzero if u and u e are 
to be nonconstant solutions, the boundary conditions give us the following conditions on a: 



2 sin(a) + a cos(a) cos(a) + a sin(a) 

— sin(a) cos(a) 

2 cos(a) — a sin(a) sin(a) — a cos(a) 
cos(a) — sin(a) 



for u a , 
for u e . 



These simplify to 

3 cos(a) sin(a) + a = for u D , 

3 cos(a) sin(a) — a = for u e . 

These can be rewritten, using the double-angle formula, as 

sin(2a) + 2a/3 = for u Q , 

sin(2a) — 2a/3 = for u e . 

It is easy to see graphically that the latter has only the solution a = 0. The former has two 
nonnegative solutions, a = and a w 1.13943. We therefore have no odd solutions in this case, 
and exactly one even solution with a w 1.13943. In particular, this means that the curves given 
by the eigenvalues plotted against r may cross the curve co = — t 2 /4 for t < only once, at 
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t re -2.5966. □ 
The hyperbolic regime: — r 2 /4 > w 

In this case, our characteristic equation r 4 — rr — u = has non-real, non-purely-imaginary 
complex roots; the real part of each of these roots then gives us a hyperbolic function in the 
solution. This is the most difficult case, and we are only able to obtain a partial result. 

Proof for the hyperbolic regime. Elementary differential equation theory gives us that the gen- 
eral solutions will be of the forms cos(ax) sinh(fex), sin(ax) sinh(6x), cos(ax) cosh(6x), and 
sin(aa;) cosh(6x), with u = —(a 2 + b 2 ) 2 and r = 2b 2 — 2a 2 ; a and b are taken to be positive as 
usual. As before we need only consider the even and odd solutions 

u (x) = Acos(ax) sinh(&ir) -I- B sin(ax) cosh(6a;) and 
u e {x) = C sin(aa;) sinh(6.x) + D cos(ax) cosh(6a;). 

Their associated boundary conditions are 

A (rcos(et) sinh(6) — 4a6sin(a) cosh(fc)) = —B (r sin(a) cosh(&) + 4a6cos(a) sinh(fe)) (7.10) 

A(acos(a) cosh(fr) + &sin(a) sinh(&)) = — _B(asin(a) sinh(fe) — 6cos(a) cosh(6)) (7-11) 

c{t sin(a) sinh(6) + 4a6cos(a) cosh(fe)^ = — D (r cos(a) cosh(6) — 4a&sin(a) sinh(b)) (7.12) 

C{a sin(a) cosh(6) — 6cos(a) sinh(6)) = — -D(acos(a) sinh(fe) + 6sin(a) cosh(fr)). (7-13) 

We are able to show that the constants B and D must be nonzero. 



Suppose B = but 4^0. The boundary conditions (7.111 and (7.10) can both be solved 



for tan(a); combining the resulting equations we obtain tanh 2 (6) = —^r- Recall r = 2b 2 — 2a 2 

Mi. 

i ni'A /-iaoo no oViAirfl Tim n rr m n -fi n /-l 4" o vi Vi ( J-t\ O _1_ ' 



is negative; thus we have tanh 2 (6) = 2 1 "' 

Set D = while C ^ 0; then by the same process as above we again find tanh 2 (6) = 2 -' ' " 



The function tanh 2 (a;) is nonnegative for all x and satisfies tanh 2 (a;) < 1 for all x. Thus the 
equation tanh 2 (6) = 2+ has no solution b > for any r < 0. Therefore is impossible to take 
B = or D = while A or C are nonzero. □ 

Remark. Applying the same process in the case of A = or C — yields instead coth 2 (fr) = 
2 + T^r, which has a solution b for every r < 0. Furthermore, numerical investigations in 
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Mathematica suggest that in these cases the boundary conditions can be simultaneously satisfied 
for many values of r. 



Cascading 

Plotting the eigenvalues of the rod under compression against r, an interesting feature appears. 
Each odd eigenvalue branch in the trigonometric regime stays "close" to a line t/j^, where /i& is a 
Neumann eigenvalue with an odd eigenfunction, and then sharply drops down and travels along 
the line r/i^ +1 corresponding to the next odd Neumann eigenvalue. Similar behavior occurs with 
the even eigenvalue branches and the corresponding even Neumann eigenvalues. This is very 
much like the cascading phenomenon exhibited by Schrodinger eigenvalues discussed in |45j . 




Figure 7.2: Eigenvalues of the free rod under compression. Blue curves correspond to eigenvalues 
of even functions; purple curves indicate the eigenvalues of odd functions. The dashed curve 
is the regime boundary uj — — r 2 /4. The green and gold curves are possible odd and even 
eigenvalue branches, respectively, in the hyperbolic regime. These curves were computed under 
the assumption that all linear combination coefficients are nonzero. 
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CHAPTER 8 



Future directions 



In this chapter we examine several problems related to the free plate isoperimetric problem. One 
can generalize the plate Rayleigh quotient to plates made of material having a nonzero Poisson's 
Ratio; unfortunately, the proofs of the theorems cannot be generalized so far. Other known 
isoperimetric inequalities for membranes can be considered for the corresponding plate problems, 
namely considering a Szego-type problem for the free plate and a PPW-type problem for the 
clamped plate. Finally, it may also be worth investigating the plate isoperimetric problems in 
non-Euclidean spaces. 



Poisson's Ratio 

One generalization of the free plate problem is to account for Poisson's Ratio, a property of the 
material of the plate that describes how a rectangle of the material stretches or shrinks in one 
direction when stretched along the perpendicular direction. Our Rayleigh quotient and work 
so far all hold for a material where Poisson's Ratio is zero. Most real-world materials have 
a € [0, 1/2], although there exist some materials with negative Poisson's Ratio. 

We will assume a € (— oo, 1) for mathematical reasons made clear below. The generalized 
Rayleigh quotient is given by 

_ / fl (l - a)\ D*u\ 2 + + r\Duf dx 

and reduces to our previous quotient when a = 0. Note also that if we imposed clamped 
boundary conditions on our plate, Fact [3] would then allow us to write the Hessian term as 
(1 — cr)|Aw| 2 . Dependence on a would thus vanish from the Rayleigh quotient. Thus, when the 
plate is clamped, Poisson's ratio has no effect on the spectrum. 

Free and supported plates with nonzero Poisson's ratio have been considered previously: see 
|38j . |34j . and [35]. Payne [35] considered arbitrary Poisson's ratio but zero tension. 

7G 



We can show coercivity of our new quotient Q for a < 1 by discarding the Laplacian term 
and repeating our previous argument: 

a(u,u) + K"||w||| 2 

= (1 - a)\\D 2 u\\ 2 + <t\\Au\\ 2 + t\\Du\\ 2 + K\\u\\ 2 

>(l-a- S)\\D 2 u\\ 2 L2 + \\Du\\ 2 l2 + (k-™-5\ ||u||| 2 . 

This argument works for each value of r provided we choose S appropriately and the constant 
K large enough. Note that when a = 1, we lose the entire Hessian term and hence coercivity; 
this fact is why we only consider a < 1. 

Following our earlier derivation, we obtain the same eigenvalue equation 



AAu — tAu = luu, 



along with new natural boundary conditions on dft: 



= 



. a u 

Mu\ dn := {l-cr)-^+aAu 
Vu\ d n := (rg - (1 - a)dW 9n (P dn [(D 2 u)n]) - ^ 



= 0. 

an 



As before, n is the outward unit normal to the surface, divan the surface divergence, and Pqq 
the projection onto the tangent space of d£l. 

The generalization to nonzero a does not change the eigenvalue equation and hence the 
general form of solutions is preserved. However, the change in the Rayleigh quotient affects the 
proof of Theorem [2j which identified the fundamental mode of the ball. This in turn affects the 
proof of the isoperimetric inequality. Below I will discuss where the proof of Theorem [2] breaks 
down, first looking at the failure of the argument for the I > 1 case, and then observing that a 
modification of the argument in the I — 0, 1 case can be made for nonzero a. 

Considering eigenfunctions of the form u = R(r)Yi as in the proof of that theorem and 
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following our previous calculations, the numerator of the Rayleigh quotient N[u] becomes 



o 



+ a£ (^ (r)+ ( d - 1 r )^)-^)Vi dr 



o 



with fc = Z(Z + cZ — 2). When <r — and considering positive tension we were able to rewrite 
the numerator of the Rayleigh quotient to show was an increasing function of k for k > 1 and 
hence an increasing function of I for Z > 1. Unfortunately, when our Poisson's Ratio is nonzero, 
we can no longer complete the square to obtain increasing functions of k. There is a cross-term 
involving both k and R", but (R") 2 involves no factor of k, so we cannot complete a square 
to transform these terms into a something that is an increasing function of k for all R. At 
present I cannot see how one could rewrite these terms to have an always-increasing function of 
k; furthermore, numerical investigations with Mathematica suggest that for some choices of r 
and a, the fundamental mode for the ball corresponds to I — 2 or greater. 

The proof that the lowest eigenvalue corresponding to the eigenfunctions (ji(ar)+jbi(br))Yi 
is lower than that corresponding to (jo(ar) + r ybo(br))Yo still holds; we are able to show that 
MVi(a) has a root in (0,aoo) and that the first root of MVo(a) falls after a^. So while we are 
unable to prove the precise form of the fundamental mode, we still have established that the 
fundamental mode corresponds to some I > 1, and hence has angular dependence. Thus we 
have lost a lot of needed information about our radial function, and only know that 



P = 3l ( ar ) + lk ( br ) on [0, 1] 



for some positive I not necessarily 1 = 1. Hence our proof for the isoperimetric inequality is also 
adversely affected when a ^ 0, and I cannot see how to complete it. It is not even clear whether 
the fundamental tone should be maximal for the ball when a ^ 0. 
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Harmonic mean of low eigenvalues 

In two dimensions, Szego was able to prove a stronger statement of the Szego- Weinberger in- 
equality using conformal mappings 46, 48]. Specifically, he proved that the sum 

1 1 

1 

Mi M2 

is minimal for a disk. In other words, the harmonic mean of /ii and /i2 is maximal for the disk. 
Our investigation in Chapter [3] with the moment of inertia suggests a similar result for the free 
plate, since the moment of inertia is minimal for a ball. That is, for the free plate, we conjecture 

- V — > 1 

Curved spaces 

We have taken our region Q, to be in Euclidian space M. d , but we could consider the same 
eigenvalue problem on a region in spaces of constant curvature: the sphere and hyperbolic 
space. Other eigenvalue inequalities have been proven in these spaces [3]. In particular, the 
Szego- Weinberger inequality for was proved for domains on the unit sphere by Ashbaugh and 
Benguria [7] . Another direction of generalization would be Hersch-type bounds for metrics on 
the whole sphere or torus; see [26] . 
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APPENDIX 



Calculus facts 



This appendix collects some calculus facts used in the proof of Theorem [I] Below, p(r) is any 
smooth function of the radial coordinate and the xi are the usual rectangular coordinates. We 
also define functions 

p(r) 



u k = x k - 

r 



Fact .1. We have the 



Ew 2 = p 2 

fc=i 



fe=i 



jr\DW\ 2 = (p f r + 3{ ^(p-rpr 
fe=i 

E(^) 2 = ((«-i)^-p") 2 . 

fc=l v / 

Proof. The first is immediate from x\ + ■ ■ ■ + x 2 d = r 2 . Now note that 

dr _ x^k 
dx k r 



so that 



du k p XiXk , n 



Thus 

d 2 u k x k S l:j + XiS jk + XjS ik , XiXjX k f 3(p - rp') 



dxidxj r 3 P P + 



and 



d 2 u k x k + 2xiS ik , x 2 x k (Z(p-rp') „ 

= ^ iP-rp) + - W [—^—+P 
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Write A = (p — rp') and B = 3r 2 (p — rp') + p". Then the sum of the gradients becomes 



dd 2 



fe=l fe=l i=l 



-,2 



-(^-D^ + (P') 2 



by using the <5jfc to simplify. The sum of the Hessian terms becomes 



d d d d / r. \ z ] X \ 2 

|£) 2 U | 2 — ^""^ ^""^ ^""^ / _ x k°ij + x iOj k + x jOik _j_ x i x j x k ^ j 

fc=l fc=l i=l j=l V r r J 

= (p") 2 + 3(d-l)^, 



after collapsing the sums and simplifying. Similarly, the sum of the Laplacians becomes 
V^a ^2 V* f 1 + 2S ik x\x k ^ ( 1 + 25 k . 3 x)x k \ 

k=l i,j,k=l V / \ / 

The next two facts express the norm of the Hessian matrix in terms of the Laplacian and a 
divergence. 

Fact .2. Ifu(x) is C 3 -smooth, then 

\D 2 u\ 2 = - (A\Du\ 2 - D{Au) ■ Du - D{Au) ■ Duj 
= X -D ■ (p\Du\ 2 - AuDu - AuDu^j + | Au\ 2 . 

Here u could be complex-valued, with u denoting its complex conjugate. 
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Proof. We have 



id 2 



D ■ (p\Du\ 2 - AuDu - AuDu^j + 2\Ai 

= A\Du\ 2 - D(Au) ■ Du - AuAu - D(Au) ■ Du - AuAu + 2\Au\ 2 

= A\Du\ 2 - D(Au) ■ Du - D(Au) ■ Du 

d 

= y ] (%u XkXl u XkXl + u Xk u XkXlXl + u XkXlXl u Xk ^j 

k,l=l 

d d 

~ y ] u x k x k xiU X i — ^ ] u x k x k xi u xi 
j,k=l J,k=l 

= 2\D 2 u\ 2 . □ 



Fact .3. Ifue C 3 (n) and Du = on dQ, then 

D 2 u\ 2 dx= / \Au\ 2 dx. 
n Ja 

Proof. Integrate Fact [2] and apply the divergence theorem. □ 
Thus if u satisfies clamped boundary conditions u — du/dn — 0, then the Hessian term in 



the Rayleigh quotient (1.1 1 can be replaced with |Au 



2 



t s2 
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